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1.  INTRODUCTION 


Although  knowledge  of  the  amount  of  solar  radiation  absorbed  and  re¬ 
flected  back  to  space  by  the  ground  and  atmosphere  is  pertinent  to  a  large 
variety  of  geophysical  processes,  the  application  for  which  the  present 
procedure  has  been  formulated  involves  the  heat  balance  of  terrestrial  ob¬ 
jects  and  backgrounds.  The  application  requires  a  simple  means  to  approx¬ 
imate  the  flux  of  both  upward-  and  downward-directed  solar  radiation  through 
the  atmosphere  using  only  routine  surface  weather  observations.  These  re¬ 
quirements  have  served  to  limit  and  focus  the  structure  of  the  procedure. 
Although  the  downward  flux  of  solar  radiation  at  the  ground  is  the  only 
parameter  which  is  specifically  required,  this  flux,  as  well  as  each  of 
the  fluxes  at  all  levels,  is  dependent  upon  all  of  the  fluxes.  This  means 
that  it  is  necessary,  at  least  implicitly,  to  know  or  to  solve  for  all  of 
the  fluxes  in  order  to  obtain  any  of  the  fluxes.  Thus,  in  obtaining  the 
downward  flux  of  solar  radiation  at  the  ground,  one  also  has  available  the 
amount  of  radiation  absorbed  in  the  atmosphere  and  reflected  back  t  ; 
space . 

The  author  had  developed  (Shapiro,  1972)^  a  simple  model  for  the  flux 
of  solar  radiation  through  the  atmosphere  wh.ich  was  readily  adaptable  for 
the  present  purpose.  The  model,  however,  requires  knowledge  of  the  coef¬ 
ficients  of  reflectivity,  transmissivity  and  absorptivity  of  each  of  the 
n,  plane-parallel  layers  into  which  the  atmosphere  is  divided.  Further¬ 
more,  since  these  coefficients  are  most  sensitive  to  variations  in  cloud 
type  and  amount,  knowledge  of  the  absorption,  reflection  and  transmission 
of  solar  radiation  for  individual  cloud  types  is  particularly  necessary. 
Ultimately,  such  knowledge  must  be  based  up>on  direct  observation.  Perhaps 
the  most  satisfactory  approach  is  the  direct  observation  of  reflection, 
transmission  and  absorption  of  individual  clouds  from  instrumented  air¬ 
craft.  Unfortunately,  such  measurements  are  difficult  and  expensive  and 
there  exists  today  only  a  limited  number  of  such  measurements  for  a  lim¬ 
ited  number  of  cloud  types.  Although  there  have  been  some  recent  measure¬ 
ments  (Paltridge  and  Platt,  1981";  Reynolds  et  al.,  1975^,-  Paltridge,  1973*^ 
and  1974  ;  Cox  et  al.,  1973  Drummond  and  Hickey,  1971  )  ,  and  some  not  so 

Because  of  the  number  of  references  to  be  included  as  footnotes  on  tins 
page,  the  rfiaac"  is  referred  to  the  list  of  references,  page  5u. 
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recent  (Feigel'son,  1966  Predoehl  and  Spano,  1965';  Hanson  and  Viebrock, 
1964^^) ,  much  of  the  available  information  was  obtained  more  than  20  years 
ago  (Roach,  1961^^;  Robinson,  1958^^;  Fritz  and  MacDonald,  1951  ;  Fritz,- 

1950^"^;  Neiburger,  1949^^)  . 

In  addition  to  these  measurements,  there  have  been  some  direct  measure- 

IS 

ments  of  cloud  microphysics,  such  as  those  of  Reynolds  et  al.  (1978)  , 

Platt  (1976)^^,  Paltridge  (1974)^,  among  others.  These  measurements  of 

drop-size  distribution  and  liquid  water  content  have  served  as  the  basis  for 

a  large  number  of  theoretical  computations.  Much  of  the  present  knowledge 

of  the  flux  of  solar  radiation  in  the  atmosphere  and  its  variation  under  a 

variety  of  conditions  is  due  to  such  computations,  which  are  based  both  upon 

scattering  theory  and  Monte  Carlo  calculations.  The  results  of  course  are 

dependent  upon  the  modeling  assumptions  which  are  made,  but  these  in  recent 

years  have  become  increasingly  sophisticated  and  presumably  increasingly 

realistic.  Prominent  among  these  studies  is  the  recent  monograph  by  Welch 
18 

et  al.  (1980)  which  contains  a  large  number  and  variety  of  computations. 

19 

Other  recent  reviews  are  those  of  Van  De  Hulst  (1980)  ,  Fouquart  ot  al. 

(1980)''^  and  Lenoble  (1977)^^.  An  indication  of  the  volume  and  scope  of 

such  studies  may  be  given  by  a  sampling  of  some  recent  contributions: 

'’2  23 

Newinger  and  Bahnke  (1581)  ;  Platt  (1981)  ;  Welch  and  Zdunkowski  (1981a, 

24  yf>  27  2P 

b)  '  j  Fouquart  and  Bonnel  (1980)  ;  Leighton  (1980)  ;  Manton  (1900)  j 

2n  30 

Meador  and  Weaver  (1980)  ;  Zdunkowski  ot  ol.  (1980)  ;  Davis  ct  al.  (1979a, 

313?  33  34 

b)  '  ‘ ;  Liou  and  Wittman  (1979)  ' ;  Schaller  (1979)  ;  Schmetz  and  Raschke 

3S  3637  38 

(1979)  ;  Stephens  (1976,1979)  '  ;  Davies  (1978)  ;  Kerschgens  et  al. 

39  40  41  4"’ 

(1978)  ;  Twomey  (1976,1977)  ;  Wendling  (1977)  Wiscombe  (1976a, b, 

43  44. 4S  46  47  48.49 

1977)  '  ;  Liou  (1974,1976)  '  ;  McKee  and  Cox  (1974,1976)  ;  Welch 

ot  al  (1976)^^;  Wiscombe  and  Grams  (1976)  Kdunkow.ski  and  Korb  (1974)^^. 

Although  there  .ire  r.izeable  differences  among  the  measuremenrs  and 

among  the  calculations  for  similar  cloud  conditions,  the  results,  insofar  as 

comparisons  are  feasible,  are  generally  individually  and  mutually  consistent. 

The  central  feature  cif  both  the  studies  .md  the  direct  measurements  is  the 

not  unexpected  result  that  the  fractions  of  sunlight  reflected,  transmitted 

and  absorbeci  are  proportional  to  the  density  and  thickness  of  the  clouds, 

or  in  other  words,  are  functions  of  cloud  type. 

53 

llaurwitz  (1948)  ,  by  means  of  hourly  cloud  observations  and  coordin¬ 


ated  solar  insolation  measurements,  was  able  to  obtain  estimates  of  mean 


transmissivity  as  a  function  of  cloud  type  and  solar  zenith  angle.  Although 

there  have  been  a  number  of  similar  studies  in  more  recent  years,  such  ar 

Kasten  and  Czeplak  (1980)^^,  Atwater  and  Ball  Tabata  (1964)^^, 

57  58 

Lumb  (1964)  ,  Vowinckel  and  Orvig  (1962)  ,  Haurwitz's  extensive  study  re¬ 

mains  definitive  for  its  type.  Nevertheless,  Haurwitz's  results  were  ob¬ 
tained  from  only  eight  years  of  data  at  Blue  Hill,  Massachusetts  and  are  defi- 
cient  in  several  cloud  types,  particularly  for  the  larger  zenith  angles. 

'I’heorct  ica i  model  calculations  have  the  adv.intage  of  completeness  in 
that  results  can  be  obtained  for  any  assumed  set  of  cloud  conditions.  How¬ 
ever  those  models  ^^re  not  suitable  for  routine  computations  since  they  re-- 
quire  cloud  microphysical  information  which  is  not  routinely  available. 
Furthermore,  there  have  been  too  few  coordinated  measurements  of  cloud  type 
and  cloud  microphysics  to  substitute  cloud  type  for  microphysics  in  these 
models.  The  present  approach,  which  has  been  deliberately  kept  simple,  is 
desiqtiod  for  appl  ii:a  t  i  on  witli  information  available  solely  from  routine  sur¬ 
face  weather  observations.  It  makes  use  of  a  very  simple  two-stream  ajiprox- 
Ltn.ition  and  moan  climatological  reflection  and  transmission  coefficients 
determined  from  the  HOIMFT  (1977?^  data  tabulations. 

Atwater  and  Ball  (1978,1981)^^'^’^  have  also  proposed  a  simple  model  for 
the  com)iutation  of  radiation  received  at  the  ground  based  upon  standard  sur¬ 
face  meteorological  ooservations .  It  is  believed  that  the 

model  incorporates  multiple  reflections  and  the  observed  nonlinoarities  of 
reflectance,  transmi ttance  and  absorptance  with  fractional  cloud  amounts 
in  a  more  realistic  fashion. 

S4 

Kasten,  F.  and  G.  Czeplak,  1980:  Solar  and  terrestrial  radiation  dependent 
on  the  amount  and  type  of  cloud.  Solar  Energy,  24,  177-190. 

65 

Atwater,  M.  A.  and  J.  T.  Ball,  1978:  A  numerical  solar  radiation  model  based 
on  standard  meteorological  observations.  Solar  Energy,  21,  163-170. 

S(i 

Tabata,  S.,  1964:  Insolation  in  relation  to  cloud  amount  and  sun's  alti¬ 
tude.  Studies  in  Oceanography,  Hidaka  Volume,  202-210. 

57 

Lumb,  F.  E.,  1964:  The  influence  of  cloud  on  hourly  amounts  of  total 

solar  radiation  at  the  sea  surface.  Quart.  J.  Roy.  Meteor.  Soc . ,  90, 
43-56. 

58 

Vowinckel,  E.  and  S-  Orvig,  1962:  Relations  between  solar  radiation  and 
cloud  type  in  the  Arctic.  J.  Appl .  Meteor.,  1,  552-559. 

59 
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2.  MODEL  STRUCTURE 


The  model  atmosphere  is  composed  of  n,  plane-parallel  homogeneous  lay¬ 
ers  of  arbitrary  thickness.  In  addition,  the  ground  surface  constitutes 
another  "layer"  with  the  special  property  that  its  transmission  is  zero. 

Of  the  radiation  incident  on  layer  k,  the  fractions  reflected,  trans¬ 


mitted  and  absorbed  are  given  by  R^^,  Tj_,  and  respectively,  where 


R  +  T,  +  A,  ==  1. 
k  k  k 


(1) 


For  simplicity,  R  ,  T  and  A,  are  assumed  not  to  depend  upon  the  di- 

K  K  iC 

rection  of  the  incident  flux;  that  is,  they  are  assumed  to  be  the  same  for 
both  upward  and  downward  fluxes.  However,  as  will  be  seen,  some  distinc¬ 
tion  is  made  between  direct  and  diffuse  radiation,  and  to  the  extent  that 
the  direct  and  diffuse  radiative  components  depend  upon  the  direction  of 
the  flux,  so  do  the  coef f icients. 

Ai\(il  her  simplification  implicit  in  Eq.  (1)  is  that  solar  radiation  may 
bo  treated  as  quasi-monochromatic .  Selective  absorption  and  scattering  by 
uase.s  .uai  aerosols  are  troatCi.i  in  a  crude  fashion  by'  the  choice  of  values 
assigned  to  R^,  .T^^,  and  A^^  for  the  various  atmospheric  layers. 

In  Figure  I,  the  horizontal  linos  numbered  1  though  n  represent  n 
homogeneous  atmospheric  layers  and  g  represents-  the  ground  surface.  X 

)v 

represents  the  radiation  impinging  upon  layer  k+1  from  above  and  repre¬ 
sents  the  upward-directed  radiation  emanating  from  layer  k+1.  Thus  is 
the  amount  of  radiation  from  the  sun  reaching  the  top  of  the  atmosphere 
vertically  incident  on  a  unit  horizontal  surface,  amount  of  sol¬ 

ar  radiation  reflected  back  to  space  by  t^o  oarth-atmor.plicre  system.  X^, 
the  quantity  in  which  we  are  particularly  interested.,  is  the  solar  radia¬ 
tion,  both  direct  and  diffuse,  received  at  the  ground.  Each  X  and  Y  is 

K  K 

composed  of  differing  fractions  of  direct  and  diffuse  radiation,  and  while 
it  would  be  possible  to  take  account  of  the  direct  and  diffuse  radiation 
j;ep.irato  ly ,  it  is  lelt  that  this  would  needlessly  comj.-'l  icate  the  procedure, 
v'onseepient  1  y ,  a  simple  but  arbitrary  distinction  is  made',  whif:h  wiiile 
crude,  should  account,  at  least  in  part,  for  the  difference  between  the 
direct  and  diffuse  components.  All  radiation  is  considered  direct  except 
for  radiation  transmitted  through  a  layer  containing  at  least  7/8  of  any 
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Figure  1.  Flux  of  solar  radiation  through  an  atmosphere 

consisting  of  n  homogeneous  layers  and  a  ground 
surface . 


thick  cloud.  All  radiation,  once  transmitted  through  such  a  cloud  layer, 
is  considered  diffuse.  All  clouds,  except  for  thin  cirrus  and  cir rostratus , 
are  considered  thick. 

The  aibove  system  is  described  by  2n+2  linear  equations  for  the  2n+2 

fluxes  ...,  X  ;  Y_,Y  ,  ...,  Y  )  in  terms  of  the  specified  coeffi- 

0  1  n  u  J  fi 

cients  K  and  T,  ,  (k  =  0,  1,  . . . ,  n) .  The  sequence  of  equations  has  the 
K  K 

form 

\  \  \  'k 


---  R. 


k+1  k 


+  T 


k+1  k+1 


(3) 


where,  of  course,  =  1 ,  =  0 ,  R  , 

'  0  0  I)  + 1 


R  ,  and  T  =  T  =  0. 
q  n+ 1  q 


Eq.  (2)  states  that  the  downward  flux  of  r.idiation  leaving  any  layer 
k  is  equal  to  the  fractional  transmission  of  that  layer  (T^^)  times  the 
downward  flux  of  radiation  reaching  that  layer  from  above  (X^  plus  the 
fractional  reflection  of  that  layer  (R^,)  times  the  upward  flux  of  radiation 
teaching  that  layer  from  below  (f^.)  .  Eq.  (3)  makes  comiiarable  statements 
about  the  upward-directed  flux. 

The  absorption  of  solar  radiation  in  any  layer  k  is  obtained  from  Eqs. 
(2)  and  (3)  and  is  given  by 


Vi  ^ 


'“’k 


k-1 


(4) 


and  the  total  amount  of  energy  absorbed  in  the  atmosphere  is 

S  =  f  S  =  ^0  -  "O  -  "n-  <=> 

The  system  of  Eqs.  (2)  and  (3)  can  be  solved  lor  any  of  the  t luxes 

in  terms  of  specified  values  for  the  reflection  and  transmission  coeffi¬ 
cients.  In  general,  such  solutions  require  the  inversion  of  a  (?n+2)  square 
matrix,  a  not  inconsiderable  task  for  even  moderately  larqe  values  of  n. 
Fortunately,  because  tne  matrix  contains  only  a  limited  nui.iber  of  non-zero 
elements,  it  was  possible  to  find  closed  form  solutions  for  Y^  and  for 

any  n.  Ihis  imi'lies  closed  form  solutions  for  all  of  the  fluxes,  since 
they  can  be  obtained  directly  by  a  simple  step.wise  j.rocuss  from  the  sys¬ 
tem  of  Equ.  (21  and  (3)  once  a  solution  for  or  X_^  is  known.  The  solu- 
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tions 


wliere 


The  particular  solutions  of  interest  in  the  present  application  are 
for  n  =  3.  This  is  the  simplest  geometry  that  makes  use  of  the  standard 
cloud  code  information  which  categorizes  clouds  into  high,  middle,  and  low 
cloud  types.  These  solutions  are; 

''o  -  [“l  *  Vl\''  *  "3  "dj'"  \ 


X3  .  TjTjTjX^/D^ 


w!'.ere 


h  "  -  h“2>"  -  "2"3>  - 

°2  “  '>3'‘'l‘^2  -  “l''3'"2"'  -  ‘‘l'‘2V3"  ^  ‘Vs'" 


+  +  B-R-^R  +  R,R.,R.,R  +  R,R,^R  +  R^R^^^R  T,^) 

123  Z.tq  123g  13c2  ]?g3 


-  R.R.'^R^^R  . 

1  2  3  q 


It  is  interesting  to  note  that  the  solution  for  is  expressed  as  a 
sum  of  terms,  each  of  which  contains  the  contribution  from  a  separate  lay¬ 
er.  The  interactions  among  the  various  layers  for  both  Y  and  X  are  con- 

0  n 

tained  in  the  expressions  for  and  . 

As  an  aid  in  the  visualization  of  the  moaning  of  the  solutions.  Fig¬ 
ure  2  traces  each  individual  component  of  the  radiative  fluxes  for  a  3- 
Inyer  atmospliere  in  which  the  ground  surface  reflectivity  r  =  o  In  this 

Q 

cas(x  i:qs.  (8;  ,iiid  !0)  arc  r.  imp-l  i  f  ied ,  since  in  addirinn  to  rc^movinu  the 
term  ccmtaininq  R  in  Kg.  (H),  U,,  in  Kqs.  (B)  and  ('))  rcducrr-s  to  P  . 

U  .i  1 

figure  2  starts  with  the  ext ra-terrestr i al  radiation  X^  impinging  on 
the  top  of  layer  1.  Part  of  this  radiation  (Rj^q)  is  immediately  reflected 
back  to  space  and  part  (T^X^)  is  transmitted  to  layer  2.  I'art  of  the  lat¬ 
ter  component  is  transmitted  to  layer  3  part  reflected  up  to 
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the  bottom  of  layer  1  ‘  this  latter  component  (^2’^!  ^0^ 

transmitted  through  layer  1  to  space  and  part  (R^R^T^^X^)  is  returned  to  the 
top  of  layer  2.  This  returned  component  is  ^^^^2  component 

which  had  previously  been  transmitteu  to  the  top  of  layer  2.  Therefore,  the 
next  component  transmitted  to  space  is  ^,^2  antecedent  component; 

namely,  ^1  ^0*  next  comijonent  transmitted  to  space  arising  from 

the  earlier  component  and  each  fol’ owing  comi>onent  is  times  its 

antecedent  component.  This  infinite  series  is  a  geometric  progression, 

Vl\  ^  •  •  •  -  ''2"’l^^0  /  ‘^1- 

Similarly,  an  infinite  series  of  flux  terms  emanating  from  I^^X^  is  trans¬ 
mitted  down  through  layer  2  to  the  top  of  layer  3.  This  series  also  forms 
a  geometric  progression  ^  ^  j^l^2j  ^1^2^0 

[^1^2-)*^  ^l'^2^0^  whose  sum,  as  p  approaches  infinity,  is  ^  *^1' 

is  labeled  for  convenience. 

The  second  row  of  Figure  2  starts  with  .A^ ,  the  sum  of  the  radiative 
flux  emanating  from  Tj^Xq  and  reaching  layer  3  from  above,  and  traco.s  the 
Ijrogrcss  of  thi.s  radiation  as  i>ait  of  it  is  transmitted  through  layer  3  to 
the  ground,  where  it  is  absorbed  and  part  of  it  is  transmitted  up  through 
layer  2  to  the  botto, .  of  layer  1.  Those  comixinents  form  additional  infin¬ 
ite  series  whose  sums  are  /  d.,  =  C  for  the  radiation  reaching  the 

ground  and  ''  '^2 


for  the  radiation  transmitted  U))  to  layer  1, 


Again,  startin<!  witii  li^,  another  i)aii  i.if  infinite  senes  of  components 
i.s  formed.  Tile  sum  of  the  series  consisting  of  the  radiation  transmitted 
to  .space  is  the  sum  of  the  scries  consisting  of  tnc  rad¬ 

iation  transmitted  down  throucfli  layer  2  to  the  to|i  of  layer  3  is 
1' 

sum  of  the  terms  transnuttcU  lhrou.)h  layer  3  to  the  (irourid  is  "  " 

and  the  sum  ol  the  terms  tiansmitted  through  layer  2  ui<  to  the  bottom  oi 

forms  a  pair  of  infinite  sums,  ^ 


^'*-2%  /  ^'1 


O' 


From  Fq  another  p.iir  of  infinite  set  ior  is  icirmed.  The 


layer  1  is 


(out  to  sijaoe)  and  "  ‘*0  Aiycr  1).  J.,  form-. 


0 


a  pair  ol  inliiiite  series,  'ip^.I^/d^  =  K  (down  to  the  giouiui)  ,  .md  3'*'2' '  (/ ''l 
=  (up  to  the  bc'ttnm  ol  luyer  ]).  It  is  apparent  li.it  this  pioccs.-;  of 

ffirmation  of  a  pair  ol  iiil  iiiite  sums  from  e.uli  pievimis  in(  inite 
tiiiucs  indefinitely. 

If  wo  add  the  111  finite  sums  reaching  the  grouiid,  we  have 


1  •! 


sum  coii“ 


X^-C  +  H  +  K+  .  .  . 


where 

C  =  TjA^/dj  = 

H  =  T^Tp/d^  =  RiT^T^B^/d^d^  =  Rj^R3T2^T^AQ/d^d2^ 

=  RjR3T^T2\^XQ/(d^d2)^,  and 

2  2 

K  =  T3dQ/d2  =  ^  ’’3^0^'^l‘^2 

-  ^"«3"^'^/''3V'^'‘2''‘• 

C,  H,  K,  .  .  .  are  each  the  sums  of  infinite  series,  but  their  sum  (X^)  is 

an  infinite  series  of  infinite  sums,  whose  sum  is  given  by 

C  *  I.  ♦  K  *  .  .  .  •  *  ''l'3''2'^V2  •  t«i'<3>S‘'^“*l"2''  *  •  • 

12 

=  ( _ L_ _ ) 


3^0 

"A  -  ^"3^2 


2  =  ^ 


which  is  what  is  obtained  from  Eq,  (9)  when  =  0. 

Similarly,  is  an  infinite  sum  of  infinite  sums  such  that, 

'^0  =  "l"0  ^  «2^1  V'^l  ^  ^  • 

where  E  =  T^B^/d^  =  Vl'^2V’^i‘^2  '  ‘^3‘V2’^ 

1  =  T^G^d^  =  R3T^T2Fo/d^d^  =  V^l'''2 


and  where 

£  +  I  +  .  ,  .  form  the  series 
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2 


E  +  I  +  , 


’’o 


1  + 


‘'l‘'2 


‘Vj’S" 

+ - 2“  ^ 

'W 


\\ 


1  - 


'^1^2 


R3(TjT2)  Xq 


‘'l‘^2 


di(did2  -  RiR3T2')' 


Thus 


"o  ^  "o 


Vl 


R3(T^T2)' 


d^(d^d2  -  R1R3T2  ) 


which  is  what  is  obtained  from  Eq.  (8)  when  R  =  0, 


3.  REFLECTIVITY,  TRANSMISSIVITY  AND  ABSORPTIVITY  l:OEI1' U  1  ENTS 

It  is  apparent  that  the  main  impediment  to  the  inunediate  application 
of  the  model  is  a  lack  of  information  on  appropriate  values  of  the  reflec¬ 
tivity,  transmissivity  and  absorptivity  (R,  T,  A)  coefficients.  These  co¬ 
efficients  are  not  constants  but  are  functions  of  space  and  time.  We  cannot 
expect  to  account  for  all  of  the  variability,  but  much  if  not  r„ost  of  the 
variability  is  a  function  of  cloud  amount  and  type.  The  solar  zenith  angle 
is  another  important  source  of  variability.  These  sources  can  easily  be 
incorporated  from  routine  observations. 

If  the  individual  layers  are  thin  enough,  R,  T,  and  A  for  any  partic¬ 
ular  layer  can  be  treated  as  a  constant  for  that  layer,  time,  and  location. 
It  is  not  possible  to  S[)ccify  how  many  layers  the  dtmos[)here  should  be 
divided  into  in  order  to  ^u.stify  the  above  assumption,  but  it  would  un¬ 
doubtedly  be  more  than  the  three  layers  wo  will  be  using  for  our  first 
approx imat ion , 


IG 


Although  som 5  effort  has  been  made  to  measure  R,  T  and  A  directly,  as 
indicated  in  the  Introduction,  such  measurements  are  difficult  and  expen¬ 
sive,  Therefore,  too  few  cloud  types,  for  relatively  few  varieties  of  con¬ 
ditions,  have  been  measured  to  permit  confidence  in  the  results.  Further¬ 
more,  measurements  for  the  same  cloud  type  show  large  variations.  Part  of 
the  variability  is  due  to  the  natural  variability  in  cloud  microphysics  for 
the  same  clouti  type  from  time  to  time  and  from  place  tj  place.  However, 
other  sources  of  variability  almost  invariably  conta'.iinate  the  results. 
These  include  variations  in  the  reflectivity  of  the  surface  (either  ground 
or  cloud)  underlying  the  cloud  layer  being  measured  as  well  as  instrumental 
errors.  Furthermore,  measurement  of  any  (R,  T,  A)  coefficient,  such  as 
the  reflectivity  of  layer  k+1,  involve?  more  than  the  relatively 
simple  measurement  of  V.  A,  ,  the  albedo  of  la/er  (1^+1). 

K  K 

We  have  from  Eq.  (3) 


Y  T  Y 

-  _ii  k+1  k+l 

A+1  "  X,  '  X, 
k  k 


(10) 


This  means  that  the  upward  radiation  at  the  bottom  of  layer  k+1,  must 

also  be  measured  simultaneously  and  there  must  be  knowledge  of  the  trans¬ 
missivity  of  the  cloud  •  Btit  from  Eq.  (2)  we  have 


_  ^k+1  \+l^k+l 

\+l  "  X  ~  X 

k  k 


(11) 


This  means  that  one  cannot  simply  measure  ^nd  separately,  but  that 

both  must  be  measured  simultaneously,  which  means  that  four  flux  terms  Xj^ , 
y  ,  X  ,  and  y  must  be  measured  simultaneously.  However,  experiments 

JC  jC  "+  i  K  ^  i 

designed  to  measure  cloud  reflectivity  typically  measure  simultaneously  only 
the  downward  and  upward  directed  fluxes  above  the  cloud  layer.  Furthermore, 
aircraft  measurements  suffer  from  the  difficulties  inherent,  in  the  insta¬ 
bility  of  the  platform. 

If  layer  (k+1)  is  a  thick  undercast  layer  situated  above  a  relatively 
low-albedo  surface,  then  both  the  and  will  be  small  and  the  frac¬ 

tion  V  /X  will  closely  approximate  B.  , .  However,  if  the  cloud  undercast 
is  thin  or  has  thin  sf^ots  or  is  situated  above  a  relatively  higli-albedo 

surface,  then  Y  /X  may  be  a  poor  approximation  of  Furthermore,  the 

K  K  K  +  i. 
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measurements  are  likely  to  suffer  from  large  variability  under  such  circum¬ 


stances. 

Because  of  the  incomplete  and  unreliable  nature  of  the  direct  measure¬ 
ments,  we  have  obtained  the  (R,  T,  A)  coefficients  indirectly  using  the 
SOLMI’T  (1977)  data  base.  The  construction  of  the  SOLMET  data  base  was  un¬ 
dertaken  by  the  National  Climatic  Center  with  support  from  the  Department  of 
Energy.  It  involved  the  rehabilitation  of  the  historical  hourly  solar  rad¬ 
iation  data  for  26  National  Weather  Service  stations  and  the  incorporation 
of  these  data  along  with  the  standard  hourly  meteorological  surface  observa¬ 
tions  for  these  stations  into  a  standard  format  on  computer  tape.  The  raw 
hourly  global  irradiance  data  from  Eppley  bulb-type  pyranometers  contained 
calibration  errors  and  changes  in  radiation  scale  and  suffered  from  decolor- 
ization  of  the  black  coating  (Hoyt,  1979)^^.  h  goal  of  the  reliabilitation 
Was  to  provide  a  homogeneous,  consistent  set  of  data  with  at  most  a  5  percent 
error.  There  may  be  some  question  whether  this  goal  was  reached  (Hoyt,  1979) 
but  it  was  at  least  approached.  Nevertheless,  because  of  the  large  amount  of 
data,  random  errors  should  be  virtually  absent  from  mean  values.  Therefore 
it  was  felt  that  mean  values  of  global  irradiance  for  specified  sky  states 
and  zoiiith  ancles  would  be  reliable.  The  observed  irradiance  data  are  equiv¬ 
alent  to  in  the  3-layer  model. 

We  made  use  of  observed  mean  values  of  for  specified  conditions  to 
obtain  the  (R,  T,  A)  coefficients  for  each  of  the  three  layers.  In  es¬ 
sence  we  inverted  the  model  solutions  (X^)  in  order  to  obtain  the  (B,  T,  A) 
coefficients.  In  doing  this  we  have  worked  primarily  with  uniform  sky 
states;  that  is,  where  all  three  layers  are  clear  or  where  the  first  layer 
observalile  from  the  ground  is  overcast.  In  the  case  cf  clear  skies  or  high 
clouds  there  is  no  ambiguity;  however,  in  the  case  of  middle  or  low  layer 
overcast  we  have  no  knowledge  of  the  presence  or  absence  ot  higher  clouds. 
This  uncertainty  contributes  to  some  dispersion  in  the  results  for  middle 
and  low  clouds,  but  the  final  (B,  T,  A)  coefficients  which  are  obtained  re¬ 
flect  climatological  mean  conditions.  Thus,  we  do  not  claim  that  the  coef¬ 
ficients  obtained,  say,  for  some  low  cloud  type,  represent  the  characteris¬ 
tics  of  that  type,  but  rather  that  they  reflect  the  characteristics  of  that 
type  modified  by  the  climatology  of  middle  and  high  clouds  associated  with 


Hoyt,  D.  V.,  1979:  An  error  in  the  rehabilitation  of  the  National  Weath¬ 
er  Service  solar  radiation  data.  Solar  Energy,  23,  5S7-559. 
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that  type,  it  will  become  apparent  that  this  source  of  error  is  a  relatively 
iianor  consideration. 

Table  1  sunut\arizes  the  observed  data  for  the  13  stations  designated  in 
Figure  3  as  the  development  stations.  The  data  are  given  as  a  function  of 
certain  uniform  sky/weather  states  and  the  cosine  of  the  zenith  angle,  cos  z. 
The  transmission  of  solar  radiation  to  the  ground  is  expressed  as  1000  times 
the  average  fraction  of  the  extra-terrestrial  radiation  (X^) .  The  cos  z 
intervals  range  t.05  from  the  designated  central  values.  Although  nine  dif¬ 
ferent  cloud  types  are  recognized  by  the  standard  synoptic  code  at  each  of 
the  three  cloud  levels,  many  of  these  cloud  types  are  rare.  Nine  cloud 
types  occur  with  sufficient  frequency  to  obtain  reliable  mean  values.  These 
are  thin  and  thick  cirrus  (Ci)  ,  thin  and  thick  cirrostratus  (C.'s)  ,  altostratus 
(As)  ,  altocumulus  (Ac)  ,  cuni.-Tus  (Cu)  ,  stratocumulus  (Sc)  ,  and  stratus  (St)  . 

The  mean  values  of  observed  with  certain  cloud  types  were  indistinguish¬ 
able  from  each  other.  Therefore,  in  Table  1  Ci  and  Cs  are  combined  into  a 
single  high  cloud  type  (Ci/Cs) ,  As  and  Ac  are  combined  into  a  single  middle 
cloud  type  (As/Ac) ,  and  low  clouds  Sc  and  St  are  also  combined.  There  is  how¬ 
ever  a  large  difference  between  thin  and  thick  Ci.'Cs  relevant  to  tlie  observed 
mean  values  of  x^ 

In  addition  to  cloud  types  and  clear  skies.  Table  1  contains  some  weath¬ 
er  categories.  A  small  but  distinct  difference  was  found  between  values  of 
for  cases  where  fog  or  smoke  or  both  (F/k)  were  reported  and  tliose  cases 
in  which  no  present  weather  was  reported  for  the  categories  of  clear  skies 
or  high  overcast.  For  the  remaining  cloud  cateciorios  no  appreciable  differ¬ 
ence  was  discernible.  Consequently  separate  (R.  T,  A)  coefficients  were  ob¬ 
tained  for  (F,^K)  only  for  the  categories  of  clear  tnd  high  overcast.  The 
last  row  in  Table  1  contains  mean  values  of  X^  for  those  cases  in  wnicii  any 
kind  of  precipitation  is  reported.  The  moan  values  associated  with  iirecipi- 
tation  are  the  least  reliable.  This  is  partly  due  to  the  smaller  inimher  of 
cases  and  partly  due  to  tlie  fact  that  these  values  represent  weighted  means 
of  several  different  cioud  types.  Each  of  the  cbiuvi  types  may  have  some  iiat 
different  radiation  ciiaracteristics.  Tlicreforc,  witti  the  same  individual 
mean  values  for  each  cloud  tyiic,  different  distributions  of  cloud  lyjuis  will 
yield  different  observed  values  for  the  fractional  transmission. 
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SOLMET  STATION  NETWORK 


The  weighted  average  of  the  standard  deviations  associated  with  the 
values  of  Table  1  is  0.087.  while  there  is  no  consistent  variation  with 
zenith  angle,  there  is  a  pronounced  variation  with  sky  condition.  The 
average  standard  deviations  vary  by  about  a  factor  of  two  from  0.067  for 
clear  skies  to  0.132  for  overcast  low  clouds. 


TABLE  1.  OBSERVED  AVERAGE  FRACTIONAL  TRANSMISSION  TO  THE  GROUND  (X  1000) 
FOP  13  DEVELOPMENT  STATIONS.  CLOUD  VALUES  REPRESENT  OVERCAST 

_ CONDITIONS  . _ _ _ 

COS  Z 


,05 

.15 

.25 

.35 

.45 

.  55 

.65 

.75 

.85 

.95 

CLEAR 

460 

512 

619 

685 

735 

770 

783 

792 

794 

794 

CLEAR  (F/K) 

390 

427 

520 

602 

655 

6BS 

704 

715 

720 

722 

THIN  Ci/Cs 

395 

436 

526 

602 

660 

692 

712 

729 

733 

736 

THIN  Ci/Cs  (F/K) 

338 

36? 

446 

532 

591 

621 

642 

660 

666 

671 

THICK  Ci/Cs 

284 

319 

39C 

452 

500 

r.  r  r\ 

-f  ->\j 

1.*  ■' 

/  -> 

595 

r  A 

r  r\  A 

THICK  Ci/Cs  (F/K) 

257 

289 

357 

407 

456 

493 

513 

533 

540 

540 

Cu/Cb 

150 

179 

220 

260 

290 

316 

335 

3  50 

365 

380 

As/Ac 

184 

204 

244 

271 

297 

315 

325 

335 

338 

338 

Sc. 'St 

117 

1  39 

178 

208 

232 

250 

266 

280 

290 

296 

PRECIPITATION 

095 

105 

123 

1  35 

146 

154 

158 

161 

162 

162 

It  is  apparent 

that  is 

appreciably  smaller 

when 

there  is 

precipita 

tion,  especially  for 

large  values 

of  cos 

z . 

Since 

it  IS  unlikely 

that 

the 

precipitation  itself  has  such  a  large  effect  on  the  transmission  of  solar 
radiation,  it  was  assumed  for  purposes  of  calculation  that  when  precipita¬ 


tion  is  occurring  there  are  overcast  clouds  in  all  three  layers:  thick  Ci/Cs, 
As/Ac ,  and  Sc/St.  However,  when  there  is  no  precipitation  the  sky  is  as¬ 
sumed  to  be  clear  above  the  lowest  overcast  layer.  Thus  the  reflectivity 
and  absorpt ivi tv'  coefficients  for  middle  and  low  layer  overcast  are  un¬ 
doubtedly  biased  toward  higher  values. 

The  results  in  lable  1  were  obtained  separately  for  each  of  the  13 
stations  and  for  each  month  in  order  to  accommodate  possible  geographical 
or  seasonal  variations.  However,  since  no  clear  geograp'hi cal  or  seasonal 
differences  could  be  discerned,  the  separate  results  were  combined.  Pre¬ 
sumably,  if  tile  data  were  less  noisy  some  small,  but  real  geographical  and 
seasonal  differences  might  have  Ijeen  uncovered.  Tlie  data  were  also  exam¬ 
ined  separately  for  hours  before  and  after  noon,  but  here,  too,  no  system¬ 
atic  differences  were  found. 

The  values  of  Table  1  arc  shown  graphically  in  Figure  4,  where  it  is 
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possible  to  discern  more  clearly  the  variation  with  cos  z.  It  is  also  clear 
from  Figure  4  that  little  insolation  reaches  the  ground  when  precipitation 
is  occurring  and  that  this  insolation  varies  little  with  cos  z.  It  is  in¬ 
teresting  to  note,  for  cos  z  >  0.5,  that  while  the  differences  are  small, 
the  transmission  for  overcast  Cu/Cb  is  greater  than  that  for  overcast  As/Ac. 
On  the  other  nand ,  for  cos  z  0.5,  the  reverse  is  true.  A  possible  inter¬ 
pretation  of  tl'iis  result  is  that  Ou/Ob  are  convective  clouds  with  appreciable 
vertical  development.  Therefore,  even  when  the  layer  is  overcast,  there 
are  likely  to  be  small  breaks  or  thin  spots  in  the  cloud  layer.  When  the 
sun  is  high  in  the  sky  some  solar  radiation  may  teach  the  ground  relatively 
unimpeded,  but  because  of  the  vertical  cloud  development,  more  radiation 
IS  likely  to  be  blocked  by  cloud  sides  when  cos  z  <  0.5. 


(R,  T,  A)  coefficients  evaluated  for  clear  skies  and  for  the  various 
overcast  cloud  layers  are  piresenterd  in  Tables  2  and  3.  The  evaluation 
proco.s.s  was  a  stepwise  procedure  which  started  with  the  simplest  and  least 
ambiguous  sky  state  and  advanced  in  stages  to  the  most  complex  state. 

Since  it  was  our  intention  to  treat  solar  radiation  as  monochromatic  and  to 
determine  a  single  set  of  clear  Layer  coefficients,  processes  such  as  ozone 
and  water  vapor  absorption,  molecular  scattering,  and  scattering  and  absorp¬ 
tion  by  aerosols  could  be  incorporated  in  only  the  most  rudimentary  fashion. 
Since  the  largest  effect  on  insolation  is  due  to  clouds  and  since  we  are 
approximating  the  atmosphere  with  only  3  layers,  a  more  sophisticated  treat¬ 
ment  of  the  clear  layers  was  deemed  unwarranted. 

TABLF  2.  ESTIMATES  OF  HEFLEOTIVITIES ,  TRANSMISSIVITIES  AND  ABSORPTIVITIES 

_ (X  1000)  FOR  (-Di:ar  sky  layers. _ 


.05 

.15 

.25 

.  35 

.45 

.55 

.65 

.75 

.85 

.95 

D 

Rl 

104 

'^osT' 

042 

033 

028 

026 

”024 

024 

024 

82 

130 

112 

077 

056 

046 

039 
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040 

040 
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244 
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147 
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T1 

799 

823 
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927 
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937 
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7  33 

767 
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T3 
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762 

828 

861 

884 
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910 
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T3 

(F/K) 

59H 

616 

677 

74  3 

776 

794 

804 

811 

815 

818 

788 

A1 

097 

089 

074 

062 

052 

045 

042 

040 

039 

039 

A2 

137 

121 

094 

079 

065 

055 

052 

050 

049 

049 

055 

A3 

137 

121 

094 

079 

065 

055 

052 

051 

050 

050 

055 

A3 

(F/K) 

158 

153 

1  34 

110 

097 

090 
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084 

083 

082 
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TABLE  3. 


ESTIMATES  OF  REFLECTIVITIES,  TRANSMISSIVITIES  AND  ABSORPTIVITIES 
(X  lOOO)  FOR  OVERCAST  CLOUD  LAYERS. _ 


.05 


THIN  Ci/CP 

(Rl) 

240 

THICK  Ci/Cs 

(Rl) 

568 

As /Ac 

(R2) 

650 

Cu/Cb 

(R3) 

661 

Sc/St 

(R3) 

703 

THIN  Ci/Cs 

(Tl) 

663 

THICK  Ci/Cs 

(Tl) 

309 

As/Ac 

(T2) 

211 

Cu/Cb 

(T3) 

199 

Sc /St 

(T3) 

153 

THIN  Ci/Cs 

(Cl) 

000 

THICK  Ci/Cs 

(C2) 

026 

As/Ac 

(C3) 

002 

Cu/Cb 

(C3) 

003 

Sc/St 

(C3) 

007 

We  assume  that 

the 

modeled  by 


COS  Z 


.15 

.25 

.35 

.45 

.55 

.65 

.75 

.85 

.95 

D 

235 

211 

170 

143 

134 

120 

106 

103 

100 

535 

460 

401 

342 

297 

273 

249 

240 

240 

645 

621 

604 

587 

576 

569 

562 

560 

560 

560 

650 

636 

611 

592 

575 

561 

549 

535 

52U 

520 

693 

676 

660 

648 

639 

629 

620 

613 

609 

609 

676 

715 

768 

805 

817 

833 

849 

853 

857 

350 

440 

508 

578 

631 

660 

689 

700 

700 

2  30 

268 

292 

317 

333 

343 

353 

356 

356 

361 

226 

254 

287 

310 

331 

349 

363 

379 

395 

400 

173 

203 

227 

245 

258 

273 

286 

296 

302 

311 

000 

000 

000 

000 

004 

005 

005 

005 

004 

026 

026 

029 

028 

027 

025 

022 

021 

021 

004 

017 

025 

031 

036 

036 

035 

035 

035 

024 

003 

016 

023 

033 

039 

038 

037 

036 

035 

025 

013 

027 

034 

042 

048 

046 

04  3 

041 

039 

025 

reflectivity , 

transmissivity 

and  absorptivity 

are 

\ 

(12) 

^k  =  \^k 

(13) 

^k 

(14) 

where  r  ,  t  ,  and  a  are  the  corresponding  coefficients  for  clear  layers 
K  k  k 

and  D  and  T  are  the  reflectivities  and  transmissivities  for  overcast 
k  k 

cloud  layers.  c^  is  a  small  correction  added  to  the  clear  layer  absorp¬ 
tion  to  account  for  absorption  by  liquid  water  when  laver  k  is  overcast. 

r, ,  t  ,  and  a,  are  functions  of  cos  z  as  well  as  k.  0  ,  T  ,  and  c  are, 

K  K  K  K  K  K 

in  addition,  functions  of  cloud  type.  is  a  weighting  function  which 

depends  upon  fractional  cloud  cover.  A  full  discussion  of  will  be 
given,  but  for  the  present  it  is  sufficient  to  indicate  that  =  0  when 
layer  k  is  clear  and  ())|^  =  1  when  layer  k  is  overcast. 

We  assigned  tentative  values  for  the  reflectivity  stid  absorptivity 

(a^)  for  each  of  the  three  clear  layers  based  upon  the  mass  <cntained  in 


24 


each  layer,  modified  so  as  to  include  in  a  crude  fashion  ozone  and  water 
vapor  absorption  and  scatterinq  and  absorption  by  aerosols.  Care  was  taken 
to  remain  consistent  with  available  information  on  these  clear  layer  char¬ 
acteristics,  The  clear  layer  transmissivities  were  then  obtained  from 

Eq.  (1), 

H  ‘  -  'k  -  \  ‘  -  "k  -  \ 

It  can  be  seen  from  Eq.  (9)  that  depends  upon  the  ground  albedo 
in  addition  to  the  (R,  T)  coefficients  for  the  i  layers.  The  50LMET  data 
unfortunately  do  not  contain  specific  information  on  ground  albedo.  We 
therefore  used  a  climatological  mean  value  of  0.15  for  P,^  throughout  this 
study,  rather  than  attempt  to  estimate  as  a  function  of  time  and  location. 
In  practice,  however,  specific  information  on  the  state  of  the  ground  will 
permit  the  use  of  a  more  suitable  value. 

Xj  was  calculated  using  the  tentative  values  of  and  T^^  for  the 

clear  sky  layers  and  the  assumed  value  of  H  .  Relatively  minor  adjustments 

g  „ 

were  made  in  R^  and  as  necessary,  in  order  to  ensure  that  matched  the 
observed  mean  values  of  for  each  cos  z  category.  These  values  were  ap¬ 
proximately  the  same  as  those  contained  in  the  first  3  rows  of  each  part  of 
Table  2.  They  were  not  identical  because  and  for  the  atmospheric 
layers  are  not  uniquely  determined  by  X^.  Thus  it  was  necessary,  as  the 
coefficients  for  cloud  types  were  evaluated,  to  make  further  adju.stments 
in  the  coefficients  already  determined.  However,  in  all  cases  the  requited 
adjustments  were  of  a  minor  nature. 

Using  a  similar  procedure,  separate  coefficients,  R^IF/K)  and  T^(F/K), 
were  determined  for  the  bottom  layer  for  those  clear  sky  situations  in 
which  fog  or  smoke  or  fog  and  smoke  were  reported  as  obstructions  to  vis¬ 
ibility,  but  where  the  sky  is  not  obscured.  Situations  in  which  the  sky 
is  obscured  are  treated  as  a  stratus  overcast  if  the  obscuration  is  fog 
or  as  Q  precipitation  case  if  precipitation  is  the  source  of  the  obscura- 


tio.i. 

The  presence  of  smoke  or  fog  was  assumed  limited  to  the  lowest  layer. 

Therefore  in  evaluating  for  clear  skies  whi  smoke  or  fog  is  reported, 

the  same  values  of  R,  and  T.  arc  used  for  k  =  1  and  2  as  when  no  obstruc- 

k  k 

tions  to  visibility  arc  present. 

The  next  st.igc  in  the  prcjccss  was  to  evaluate  and  for  the  var- 


ious  uniform  cloud  layers.  The  general  procedure  was  to  start  with  the 
simplest  and  least  ambiguous  cases,  that  is,  with  thin  Ci/Cs  overcast,  and 
advance  to  increasingly  complex  states.  As  in  the  clear  sky  cases,  we 
used  available  data  (in  this  case  cloud  reflectivity  observations)  as  a 
first  approximation.  We  made  the  minor  additions  indicated  in  the  bottom 
section  of  Table  3  for  absorption  by  clouds  (c  )  and  determined  T,  from 
Eq.  (1)  .  In  the  process  of  matching  to  the  observed  mean  values  of 
to  obtain  the  coefficients  for,  say,  thin  Ci/Cs  overcast,  only  the  coeffi¬ 
cients  for  tlie  top  layer  (layer  1)  are  affected  by  the  presence  of  clouds, 
since  the  only  clouds  in  this  group  of  cases  are  high  thin  overcast  clouds. 
The  coefficients  for  layers  2  and  3  are  the  appropriate  clear  sky  values 
from  Table  2,  After  we  determined  the  tentative  coefficients  for  thin 
Ci/Cs,  we  tested  them  by  comparing  X^  and  X^  for  those  situations  where 
the  only  clouds  present  are  thin  Ci/Cs  and  where  smoke  and/or  fog  are  re¬ 
ported  as  obstructions  to  visibility.  These  latter  cases  were  not  used 
in  obtaining  the  tentative  coefficients  for  this  Ci/Cs,  but  since  the 
presence  of  F/K  only  affects  layer  3,  some  adjustments  in  coefficients  were 
necessary  so  that  the  observed  mean  values  of  X^  are  matched  by  calculated 
values  (X^)  whether  or  not  smoke  and/or  fog  is  present.  Inasmuch  as 
and  T^  are  appreciably  different  when  smoke  and/or  fog  are  present,  the 
requirement  for  matching  both  situations  imposes  severe  restrictions  on 
freedom  of  choice  for  and  T^  for  thin  Ci/Cs  overcast.  After  final  val¬ 
ues  for  thin  Ci/Cs  were  determined,  the  process  of  coefficient  evaluation 
was  repeated  for  thick  Ci/Cs,  including  the  requirement  for  matching  cal¬ 
culated  and  observed  values  of  insolation  for  both  the  presence  and  ab¬ 
sence  of  obstructions  to  visibility. 

The  next  stage  was  the  evaluation  of  coefficients  for  middle  layer 
overcast  (As/Ac) .  The  same  process  as  above  was  followed,  except  for  the 
part  concerned  with  (F/K) ,  since  the  observed  values  of  were  not  notice¬ 
ably  dependent  on  the  presence  or  absence  of  obstructions  to  visibility 
when  middle  or  low  overcast  was  i^resent  .  It  can  be  seen  from  Table  3  that 
the  variation  of  middle  layer  coefficients  with  cos  z  is  much  smaller  than 
that  of  the  high  layer  clouds.  This  is  obviously  due  to  the  greater  thick¬ 
ness  or  density  of  middle  clouds.  It  will  be  seen  that  the  coefficients 
for  low  clouds  also  have  a  relatively  small  variation  with  cos  z. 
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In  Tables  2  and  3,  the  coefficients  for  the  middle  and  low  layers 
have  additional  values  in  a  column  labeled  D.  This  column  refers  to  dif¬ 
fuse  radiation  coefficients  which  do  not  depend  upon  cos  z.  Distinction 
should  properly  be  made  between  direct  and  diffuse  radiation  in  all  three 
layers.  However,  to  do  so  in  a  precise  fashion  would  introduce  a  major 
complication  in  the  calculations.  Therefore,  as  indicated  previously,  we 
assume  that  all  radiative  fluxes  are  direct  except  for  radiation  below  an 
overcast  or  nearly  overcast  (cloud  fraction  >  7/B)  thick  cloud  layer,  where 
we  assume  that  all  radiative  fluxes  are  diffuse.  Thus,  constant  values 
(column  D)  are  used  for  both  layers  2  and  3  regardless  of  the  value  of 
cos  2  when  there  is  7/8  or  more  of  thick  Ci/Cs.  If  the  thick  cloud  is  sit¬ 
uated  in  layer  2,  then  diffuse  coefficients  are  used  only  for  layer  3. 

Clear-layer  D-values  (Table  2)  are  used  for  and  T^  when  there  is  a 
thick  Ci/t's  overcast  and  layer  2  is  clear.  Similarly,  in  layer  3,  clear- 
layer  D-values  are  used  when  there  is  either  a  middle  layer  overcast  or  a 
high  thick  overcast.  However,  with  regard  to  layer  3,  there  is  a  choice 
to  be  made  for  the  appropriate  clear-layer  D  values;  namely,  and  T^  or 
R^(F/K)  and  Tj(K/K).  The  latter  are  used  when  smoke  and/or  fog  are  re¬ 
ported  and  the  thick  overcast  is  a  high  cloud.  The  former  are  used  when 
no  obstructions  to  visibility  are  present,  or  regardless  of  the  presence 
or  absence  of  obscurations  to  visibility,  when  the  thick  overcast  is  a 
middle  cloud.  This  choice  is  consistent  with  the  observation,  already 
noted,  that  the  insolation  received  at  the  ground  in  the  presence  of  mid¬ 
dle  or  low  cloud  overcast  is  approximately  the  same  whether  or  not  obstruc¬ 
tions  to  visibility  are  reported. 

Overcast-layer  D-valuea  (Table  3)  are  required  for  layers  2  and  3 
when  precipitation  is  reported  since  it  is  assumed  that  under  this  circum¬ 
stance,  thick  overcast  cloud  layers  are  present  in  all  3  layers.  Further¬ 
more,  overcast -layer  D-values  are  used  in  layer  2  and/or  layer  3  when  the 
middle  and/or  low  cloud  cover  is  partial  and  is  beneath  a  thick,  higher- 
layer  overcast.  The  insolation  calculation  for  partial  cloud  cover  will 
be  discussed  below. 

The  D-values  in  Tables  2  a.id  3  were  chosen  on  the  basis  of  the  in¬ 
creased  path-length  appropriate  for  diffuse  radiation  as  compared  with  a 
unit  path-length  for  direct  radiation  (cos  z  -  1.00).  The  D-values  were 
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obtained  from  the  direct  radiation  values  by  interpolating  for  cos  z  =  0.60, 
which  yields  a  path-length  equivalent  to  1.66  times  the  unit  path-length. 

Some  minor  adjustments  were  necessary  for  the  clear  sky  D-values  in  order 
to  match  the  calculated  with  the  observed  values  of  X^.  The  overcast-layer 
D-values  were  also  adjusted  so  as  to  improve  the  agreement  between  the  com¬ 
puted  and  observed  values  of  for  the  precipitation  cases. 

The  final  stage  in  the  selection  of  (R,  T,  A)  coefficients  was  reserved 
for  the  low  cloud-layer  overcasts.  Although  the  differences  are  observed 
to  be  small,  a  proper  distinction  can  be  made  between  the  convective  clouds 
(primarily  cumulus,  out  including  some  cumulonimbus)  and  the  layered 
clouds  (stratus  and  .stratocumulus)  of  layer  2.  The  values  of  the  (R,  T,  A) 
coof f icients  in  Table  3  are  consi.stant  with  the  available  measurements  for 
those  cloud  types  and,  what  xs  probably  more  important,  yield  calculated 
values,  X^,  which  mntch  the  observations  in  Table  1  in  every  cos  z  category 
and  every  s';y  state  except  precipitation.  While  it  would  no  doubt  be  pos¬ 
sible  to  match  the  observations  with  different  sets  of  coefficients,  it  is 
unlikely  that  it  would  be  [jossible  to  do  so  with  appreciably  different  co¬ 
efficients  will'  ll  <iL  the  some  time  satisfy  the  various  internal  and  physical 
constraints  wliifh  wore  imposed.  Wo  reiterate,  however,  that  the  coefficients 
represent  tyrdcal,  temperate  latitude  average  values  in  which  a  multitude 
of  details  is  either  suppressed  or  treated  only  approximately.  Further¬ 
more,  only  LIk;  liigh  overcast  cloud  layer  coefficients  are  unambiguous. 

The  coefficients  for  middle  and  low  clouds,  when  no  precipitation  is  re¬ 
ported,  arc  evaluated  with  the  assumption  that  there  are  no  clouds  above 
tlic  lowest  overcast  layer.  This  means  that  the  reflectivity  and  absorptiv¬ 
ity  coefficients  for  the  middle  and  low  cloud  types  are  probably  somewhat 
too  high  and  the  transmissivity  coefficients  probably  somewhat  too  low. 
However,  this  bias  cannot  exceed  5  to  10  percent  since,  for  example,  the 
rui'orLed  racaeurements  for  middle  and  low  cloud  types  average  almost  pre¬ 
cisely  the  same  as  the  average  values  indicated  in  Table  3. 

The  qrc.ntost  uncertainty  with  respect  to  the  (R,  T,  A)  coefficients 
obtains  with  the  diffuse  coefficients  for  overcast  middle  and  low  clouds. 
Although  there  is  a  physical  basis,  their  choice  is  essentially  arbitrary. 
Furthermore,  the  manner  in  which  these  coefficients  are  used  in  the  model 
to  represent  the  flux  of  diffuse  radiation  is  at  best  a  crude  approxima¬ 
tion  to  the  real  atmosjihere. 
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4.  WEIGHTING  FUNCTION  FOR  FRACTIONAL  CLOUD  COVER 


In  order  to  compute  the  fraction  of  the  extra-terrestrial  radiation 
transmitted  to  the  ground  for  the  3-layer  atmosphere  (Eq.  9)  it  is  neces¬ 
sary  to  have  values  of  R^^  and  T^^  (k  =  1,2,3)  as  well  as  R^,  the  ground  al¬ 
bedo.  R^  is  specified  from  knowledge  of  the  state  of  the  ground  and  Table 
4  which  lists  values  of  ground  albedo  for  certain  frequently  occurring 


ground  conditions.  Table  4  was  compiled  from  a  number  of  different 

sources,  but  principally  from  Kondratyev  (1969,1793)^^'^^,  Robinson 

64  65 

(1966)  ,  and  the  ASHRAE  Handbook  (1977)  .  R^  and  Tj^  are  obtained 

from  Eqs.  (12)  and  (13).  In  applying  these  equations  it  is  necessary 


to  know  a  weighting  function  which  depends  upon  the  fractional 

cloud  cover  in  layer  k  as  well  as  cloud  type  and  cos  z.  Thus 


=  w  f,  (15) 

where  W  is  a  weighting  function  which  depends  upon  cloud  type,  amount 

and  cos  z.  f^  is  the  fractional  cloud  cover  of  layer  k.  As  indicated  a- 

bove,  vj,  =  0  for  f,  =  0  and  =  1  for  f.  =  1.0.  Thus  if  is  a  linear 
k  k  k  k  k 

function  of  1^^  regardless  of  cloud  type,  then  W  =  l  regardless  of  cloud 
type.  Although  this  is  usually  the  assumution  made  for  W  in  applications 
such  as  general  circulation  models  of  the  atmosphere,  it  is  not,  as  we 
shall  see,  warranted  by  the  observations,  which  show  that  0  is  not  a  lin- 

K 

ear  function  of  fj^  and  that  W  may  depart  significantly  from  unity. 

By  choosing  cloud-state  situations  in  which  only  a  single  cloud  type 
is  present,  it  is  possible  to  determine  W  from  the  SOLMET  data  for  each 
cloud  type,  cloud  fraction  and  cos  z.  We  examined  three  different  sit¬ 
uations. 


^^Kor.dratyev  k,  Va.,  1969:  Radiation  in  the  atmosphere.  Int.  Geophys. 
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,  (Ed.),  1973:  Radiation  characteristics  of  the  atmosphere  and  the 

Orth's  surface.  Russian  Trans.  NASA  TT  F-678,  .Amerind  Pub.  Co.,  New 
Delhj ,  pp .  192-223. 

^"^Robinson,  M,,  (Ed.’',  1966:  Solar  radiation.  Elsevier  Pub.  Co.,  Amster¬ 

dam,  Chap.  6,  pp .  196-221. 

^■^ASHRAH,  1977:  Handbocjk,  Fundamentals.  Am.  Soc .  Heat.,  Fefrig.  Air.  Ci.nd. 
Eng .,  N.Y.,  p.  2.9. 


TABLE  4.  ALBEDOS  FOR  VARIOUS  SURFACES 


WETNESS  UNSPECIFIED 


SOILS 

DRY 

WET 

OR  INDIFFERENT 

Dark 

0.13 

0.08 

Light 

0. 18 

0.10 

Dark-ploughed 

0.08 

0.06 

Light-ploughed 

0.16 

0.08 

Clay 

0.23 

0.16 

Sandy 

0.  25 

0.18 

Sand 

0.40 

0.20 

White  sand 

0.55 

SURFACES 

Asphalt 

0.10 

Lava 

0. 10 

Tundra 

0,20 

Steppe 

0.20 

Concrete 

0.30 

Stone 

0.30 

Desert 

0.30 

Rock 

0.  35 

0.20 

Dirt  road 

0.25 

0.18 

Clay  road 

0.30 

0.20 

VERDURE  UNSPECIFIED 

FIELDS 

GROWING 

DORMANT 

OR  INDIFFERENT 

Ta:^l  grass 

0.18 

0.  13 

0.16 

Mowed  grass 

0.26 

0.19 

0.22 

Desiduous  trees 

0.  18 

0.12 

0.15 

Coniferous 

trees 

0.  14 

0.12 

0.13 

Rice 

Beet,  wheat 

Potato 

Rye 

Cotton 

Lettuce 


0.12 
0.  18 
0.19 
0.20 
0.21 
0.  22 


SNOW 

ICE 

Fresh 

0.85 

White 

0.75 

Dense 

0.75 

Grey 

0.60 

Moist 

0.65 

Snow  and  Ice 

0.65 

Old 

0.55 

Dark  glass 

0.10 

Melting 

0.35 

WATER 


COS  Z 


.95 

.  85 

.75 

.0.5 

.  55 

.45 

.  35 

■  25 

.  15 

■  05 

DIRECT 

RADIATION 

.03 

.03 

.03 

.04 

.06 

.08 

.13 

.  23 

.41 

.  76 

TOTAL 

Calm 

.03 

.03 

.  03 

.04 

.  06 

.09 

.  1  2 

.  1  7 

.  24 

.  30 

RAIUATIUN 

Rough 

.03 

.03 

.04 

.05 

.07 

.08 

.  10 

.  12 

.  14 

.  15 

30 


(1)  When  the  fractional  cloud  is  in  layer  1,  from  Eq.  (9)  we  have 


2  2  2  2. 

D  =1-RR  -  RR  +  R»R  (R-,  ~  T  )  “  ^1  ~  ^ 

^  2  3  3g  2g3  J  ^ 


R  Xr  R  +  R,^  (T_^  -  Rt^)  +  T-5  (^-5  "  ^ 

2332  2  J12  ^ 


') 


(16) 


where  we  have  arranged  terms  so  as  to  factor  the  level  1  coefficient,  R^^. 


From  Eq.  (9)  we  also  have 


D  =  T.T  T  X  /X  =  T  T-T  /X 
2  1230o 


(17) 


vvhere  X^  is  merely  the  numerical  fraction  X^/X^^. 

Choosing  cases  in  which  there  are  no  clouds  in  layers  2  and  3,  and  in 
which  there  is  neither  precipitation  nor  obstruction  to  visibility,  we  ob¬ 
tain  the  values  of  R^,  R3,  which  are  required  in  Eqs.  (16)  and  (17) 

from  Table  2.  We  also  know  R^  and  for  both  clear  and  overcast  states. 
Eqs.  (12)  and  (13)  express  R^^  and  for  any  cloud  fraction  as  a  function 
of  the  corresponding  basic  (clear  and  overcast)  states,  in  terms  of  the  un¬ 
known  weighting  function  Eliminating  between  Eqs.  (16)  and  (17)  and 

substituting  for  R^  and  Tj  from  Eqs.  (12)  and  (13)  we  have 


Using  observed  mean  values  of  X*  from  the  13  development  stations  we  deter¬ 
mined  (D3  for  both  thin  and  thick  Ci/Cs  as  a  function  of  cloud  fraction  and 
cos  z.  Using  Eq .  (16)  we  obtained  the  weights  tabulated  in  Tables  6  and  6. 
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TABLE  5. 


VALUES  OF  THE  WEIGH'^ING  FUNCTION.  W.  FOR  THIN  Ci/Cs . 

COS  z 


r  KAui iwn 

.15 

.25 

.35 

.45 

,55 

.65 

.75 

.  1 

,375 

-.296 

-.826 

-1.405 

-.804 

.311 

.443 

-.272 

.  2 

1.669 

.516 

.  201 

.114 

.574 

.935 

.221 

-.136 

.  3 

1.334 

.344 

.  338 

.482 

.209 

.624 

.521 

.462 

.4 

1.133 

.697 

.223 

.260 

.384 

.  785 

.471 

.222 

.  5 

1.169 

.536 

.568 

-.009 

.074 

.459 

.758 

.145 

.6 

.755 

.610 

.654 

.510 

.321 

.758 

.659 

.560 

.7 

.742 

.646 

.681 

.609 

,477 

.809 

.745 

.715 

.8 

.764 

.726 

.850 

.799 

.658 

.865 

.869 

.605 

.9 

.952 

.741 

.913 

.842 

.825 

.954 

.999 

1.008 

TABLE  6.  VALUES  OF  TJ^_ WE  I^HTrNG_ FUNCTION ,  W,  FOR  THICK  Ci/Cs. 

COS  Z 

.45  .55  .65  .75 


.15 


.25 


.  35 


CLOUD  FRACTION 

.1 
^  2 
.  3 

.4 
.5 
.6 
.7 
.8 
.9  . 


1.594  1.143 
.777  .993 


.904 

.034 

.722 

.911 

.679 

.604 

.717 


.788 

.63b 

.727 

.716 

.744 

.709 

.626 


.85 


.409 

.  638 

.  326 

.778 

.  547 

.255 

.678 

.600 

.401 

.516 

.  357 

.215 

.660 

.500 

.347 

.497 

.498 

.027 

.669 

.567 

.355 

.373 

.429 

.064 

.  373 

.547 

.425 

.470 

.663 

.468 

.831 

.645 

.733 

.815 

.770 

.543 

.64  3 

.718 

.733 

.663 

.683 

.676 

.642 

.649 

.827 

.782 

.  665 

.711 

.770 

.749 

.706 

.851 

.917 

.836 

(2)  When  the  fractional  cloud  is  in  layer  2,  the  expression  for  is  a 
quadratic: 


^2  ‘*’2  ^2  ***2  ^2  "  ° 


(19) 


where 


A  =  R,  (R  -  B  r/  +  R  (P_^  -  T  ^  +  2t  t  -  20  r  +  r  "  -  t  “) 

2  13  g3  g3  2  2  22  22  2  2 


2  .  2. 


2 


2  2  ^ ^2  ”^2^ 

=  (r.  -  p.)  (R.  +  R.  -  R„R,  +  B  T,  -  R,R,R  )  +  - ;; - 

2  2  21  3g3  g3  13g  * 


2  2  2  2 
+  2R,  (R,  -  R  R,  +  R  T  )  Ip-t.  -  T-t-  +  t,  -  r.  ) 

13  g3  g3  22  22  2  2 


2  2  ‘l'^3’^2 

S  "  ^  ■  ^3%  -  ^^3  ■  ^q^3  ^  Vi  ■  ^I’^i'^g^ - ^ 


+  R,  (R,  -  R  R,^  +  R  T  (r.^  -  t.^). 
13  g3  q3  2  2 


(3)  When  the  fractional  cloud  is  in  layer  3, 
a  quadratic; 


+  B. 


+  C, 


the  expression  for  is  also 

0  (20) 


where 

^3  ^  '^q^'^2  "  '^1^2^  ^*^3^  "  ^3^  *  ^^3^3  ~  ^^3^3  *  ^3^  '  ^3^^ 


Bj  .  irj  -  Oj)  (Rj  .  R^  -  RjR/  .  RjT/  -  RjR^R^I  * 

^3 

2  2  2  2 
+  2R  (Ft^  -  R.R.  +  R,T-  )  (P.r_  -  T.t.  +  t«  -  r.  ) 

q2  12  12  3i  33  3  3 


T  T  ( t  •  T  ) 
12^3  3^ 


T  T  t 

2  2  12  3 

C  =  1  -  R,R^  --  r^(R^  +  R  -  R.R^  +  R,T.  -  R,R^K  ) - - - 

3  12  32  g  12  12  12g  * 


+  U  (R.  -  +  R,r^^)(r/  -  t^^). 

q  2  12  12  3  3 


In  both  Eqs.  (19)  and  (20)  the  quadratic  terms  are  so  small  that  one  of  the 
solutions  is  always  the  same  (to  three  decimal  places)  as  the  solutions  ob¬ 
tained  by  neglecting  the  quadratic  terms.  Therefore  there  is  no  ambiguity 
in  the  proper  choice  of  solution  for  4)^^.  Using  Eq .  (15)  we  obtained  the 

weights  tabulated  in  Tables  1,  8,  and  9. 
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TABLE  7 .  VALUES  OF  THE  WEIGHTING  FUNCTION,  W,  FOR  As/Ac. 

COS  Z 


.15 

.25 

.  35 

.45 

.55 

.65 

.75 

.85 

CLOUD  FRACTION- 

.1 

1.340 

.931 

.743 

.794 

.940 

.681 

.  696 

.471 

.2 

1.038 

.777 

.618 

.  447 

.673 

.603 

.611 

.452 

.  3 

.813 

.574 

.642 

.568 

.488 

.520 

.423 

.277 

.4 

.732 

.624 

.708 

.487 

.443 

.402 

.430 

.286 

.5 

.827 

.656 

.670 

.423 

.541 

.485 

.433 

.338 

.6 

.730 

.763 

.611 

.626 

.608 

.550 

.408 

.488 

.7 

.043 

.616 

.631 

.593 

.652 

.603 

.581 

.520 

.  8 

.718 

.782 

.644 

.630 

.  705 

.630 

.623 

.623 

.9 

.839 

.812 

.734 

.758 

.718 

.787 

.751 

.675 

TABLE  8 .  VALUES  OF  THE  WEIGHTING  FUNCTION,  W,  FOR  Sc/St . 

COS  z 


.15 

.25 

.  35 

.45 

.55 

.65 

.75 

.85 

CLOUD  FRACTION 

.  1 

2.085 

2.007 

1.450 

1.288 

1.234 

1.012 

.  716 

.  394 

.  2 

1.208 

1.140 

.  981 

.925 

.832 

.712 

.531 

.426 

.  3 

.827 

.913 

.789 

.753 

.744 

.674 

.533 

.444 

.4 

.751 

.791 

.680 

.803 

.727 

.665 

.  557 

.507 

.5 

.697 

.868 

.756 

.717 

.669 

.C26 

.579 

.504 

.6 

.628 

.800 

.626 

.636 

.679 

.649 

.585 

..522 

.7 

.817 

.782 

.653 

.712 

.633 

.647 

.640 

.554 

.8 

.763 

.857 

.666 

.7  37 

.744 

.  723 

.  640 

.653 

.9 

.861 

.744 

.  773 

.843 

.818 

.782 

.771 

.753 

TABLE  9.  VALUES  OF  THE  WEIGHTING  FUNCTION,  W,  FOR  Cu/Cb 


COS  Z 


.15 

.  25 

.35 

.45 

.55 

.65 

.75 

.85 

CLOUD  FFIACTION  — 

.1 

.533 

.690 

.  342 

.511 

.386 

.  352 

.  389 

.075 

.  2 

.717 

.564 

.408 

.307 

.258 

.359 

,  400 

.  115 

.3 

.750 

.646 

.413 

.429 

.272 

.339 

.  310 

.253 

.4 

.658 

.  565 

.489 

.458 

.384 

.371 

.  350 

.250 

.5 

.  521 

.  638 

.  368 

,401 

.465 

.  397 

.  369 

.  269 

.6 

.631 

.  579 

,  502 

.518 

.490 

.503 

.471 

.441 

.7 

.  648 

.  651 

.593 

.513 

.498 

.558 

.555 

.517 

.8 

.659 

.655 

.607 

.592 

.594 

.589 

.568 

.570 

.9 

.742 

.720 

.  692 

.643 

.655 

.652 

.676 

.630 

34 


It  is  apparent  from  Eqs,  (12)  and  (13)  that  a  weight  (W)  greater 
than  1.00  implies  that  the  effect  of  the  clouds  (in  terms  of  reflectivity 
or  transmissivity  coefficient)  is  stronger  than  would  be  indicated  by  the 
cloud  fraction  alone.  Similarly,  values  of  W  less  than  1.00  imply  that  the 
effect  of  the  clouds  is  wealter  than  would  be  indicated  bv  the  cloud  frac¬ 
tion  alone.  Values  of  W  near  1.00  imply  that  the  effect  of  clouds  is  near¬ 
ly  linear  with  cloud  fraction.  Negative  values  of  W  mean  that  more  insola¬ 
tion  reaches  the  ground  with  some  clouds  present  than  when  tfie  s)cy  is  com¬ 
pletely  clear.  There  are  a  few  negative  weights  for  small  cloud  fractions 
of  cirrus  or  cirrostratus .  On  the  other  hand,  for  this  same  cloud  type 
there  are  some  weights  greater  than  1.00,  expecially  when  the  sun  is  low  in 
the  sky.  Negative  weights  do  not  occur  with  any  of  the  other  cloud  ty['es; 
however,  cumulus  clouds  have  the  lowest  average  weight  (0.496) ,  which  is 
probably  due  to  the  fact  that  t)iere  is  frequent  reflection  of  sunlight 
down  to  the  ground  from  t)ie  sides  of  these  clouds.  Stratiform  low  clouds 
have  the  highest  average  weight  (0.781),  which  is,  nevertheless,  consider¬ 
ably  less  tlian  1.00.  Thus  it  is  apparent  that  it  is  improper  to  assume 
that  is  a  linear  function  of  cloud  fraction. 

5.  appliqation  of  the  model 

Routine  synoptic  observations  contain  all  of  the  infcirmation  needed 
to  apply  the  solar  insolation  model;  namely,  type  and  amount  of  clouds,  ob¬ 
structions  to  visibility  and  precipitation.  The  proper  choice  of  ground 
albedo  may  also  be  considered  a  function  of  [-resent  and  piast  synofitic  data 
insofar  as  they  contain  information  on  such  considerations  as  snow  on  the 
ground  or  the  wetness  of  tlie  qround.  Ttie  additional  information  needed 
for  the  choice  of  ground  albedo  from  Table  4  is  obtained  from  knowledge  of 
the  physiography  and  climatology  of  the  location. 

In  addition  to  R  ,  solution  of  Eq.  (9)  requires  knowledge  of  R, 

g  ^ 

and  T|^  for  the  three  layers  (k  -  1,  2,  3).  The  model  recognizes  nine 
basic  or  uniform  states,  three  in  layer  1,  two  in  layer  2,  and  four  in 
layer  3.  Each  layer  has  a  basic  state  in  which  the  layer  is  clear.  In 
addition,  the  bottom  layer  (layer  3)  has  a  basic  state  consisting  of  smoke 
and/or  foq  occurring  in  conjunction  with  an  otherwise  clear  layer  3.  Add¬ 
ed  to  these  four  clear-layer  basic  states,  there  are  five  overcast-layer 
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basic  states;  thin  Ci/Cs  or  thick  Ci/Cs  in  layer  1,  As/Ac  in  layer  2  and 
either  Cu/Cb  or  sc/St  in  layer  3.  All  possible  sky  states  including  frac¬ 
tional  cloudiness  are  represented  by  combinations  of  the  basic  states. 
Precipitation  is  assumed  to  be  represented  by  overcast  in  all  three  layers: 
thick  Ci/Cs  in  layer  1,  As/Ac  in  layer  2,  and  Sc/St  in  layer  3. 

The  first  step  in  the  application  of  the  model  is  to  use  the  synoptic 
information  to  obtain  and  the  appropriate  basic  state  values  of  and 
Tj^.  If  in  any  layer  k,  the  cloud  fraction  is  either  0  or  1,  then  only  one 
basic  state  value  is  needed  (r^^  or  for  R^  and  (tj^  or  Tj^)  for  T  .  If, 
however,  there  are  fractional  clouds  in  layer  k  equal  to  or  qreater  than 
0.1  but  less  than  or  equal  to  0.9,  then  both  the  clear  and  overcast  basic 
state  values  are  needed  in  order  to  determine  Che  appropriate  R^  and  Tj^. 

The  basic  state  values  for  R,  and  T,  in  Tables  2  and  3  are  tabulated 

k  K 

as  mean  values  for  each  cos  z  interval ,  which  means  that  they  are  appro¬ 
priate  only  for  the  midpoint  values  in  each  interval.  Inasmuch  as  cos  z 
will  be  observed  as  a  continuous  variable,  it  is  necessary  to  have  a  pro¬ 
cedure  for  obtaininq  values  for  the  basic  state  coefficients  for  inter¬ 
mediate  values  of  cos  z.  It  was  found  that  more  than  98  percent  of  the 
variance  of  each  basic  state  coefficient  could  be  explained  with  a  cubic 
polynomial  in  co?  z.  Therefore,  it  was  decidea  r<<  use  the  cubic  polyno¬ 
mial  to  obtain  the  basic  state  coefficients  ratner  than  storing  the  values 
of  Tables  2  and  3  in  the  computer  and  interpolating  for  intei-med iate  val¬ 
ues  of  cos  1.  The  polynomial  coefficients  for  each  basic  state  for  Rj^  and 

T,  are  listed  in  Tables  10  and  11  respectively, 
k 

To  obtain  the  basic  state  values  of  R.  and  T,  from  Tables  10  and  11 

k  k 

it  is  necessary  to  obtain  cos  z. 

cos  z  =  sin  0  sin  D  +  cos  0  cos  D  cos  h  (21) 

where  0  is  the  latitude,  U  is  the  declination,  and  h  is  the  hour  angle  of 
the  sun.  h  is  the  angle  bet.wei?n  toe  local  meridian  and  the  meridian  of 
the  sun.  It  is  zero  at  local  solar  noon  and  -15  degrees  or  +15  degrees 
foi  every  hour  before  or  after  local  solar  noon. 

Step  1  in  the  application  of  the  mcxJel  is  tcj  determine  the  amounts 
and  types  of  clouds  in  each  of  the  three  layers  from  the  synoptic  data 
and  then  to  determine  the  values  of  and  tj_  and/or  and  trom  the 
appropriate  sets  of  polynomial  coefficients.  If  any  of  the  three  layer.' 


lb 


i 


i 
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are  not  in  a  basic  state  condition  it  is  necessary  to  use  Eqs.  (12)  and 
(13)  to  obtain  and  for  that  layer.  This  means  that  we  must  deter¬ 

mine  'J>j^  from  Eq.  (15).  Again,  it  was  found  expedient  to  express  K  as  a 
polynomial  rather  than  apply  two-dimensional  interpolation  to  the  values 
tabulated  in  Tables  5  through  9.  It  can  be  seen  that  these  tables  do  not 
list  values  of  W  for  cos  z  =■  0.05  and  0.95.  In  general  there  were  too  few 
cases  for  cos  z  ^  .90  to  obtain  reliable  values  for  this  category.  Values 
of  W  for  cos  z  =  0.05  are  also  unreliable  because  of  the  small  values  of 

If 

X  in  the  denominator  of  X  which  enter  into  the  solution  of  .  By  eval- 

(J  j  K 

uatinq  the  polynomial  coefficients  from  the  more  reliable  values  of  W  but 
applying  the  polynomials  thus  determined  to  all  values  of  cos  z,  we  auto¬ 
matically  avoid  the  problem  of  unreliability  for  extreme  values  of  cos  z. 
Since  W  is  a  function  of  bot)i  cloud  fraction  as  well  as  cos  z,  we  esti¬ 
mated  the  weiglits  by  means  of  a  separate  bi-quadratic  polynomial  for  each 
cloud  type.  The  coefficients  are  listed  in  Table  12.  Tlie  fraction  of  the 
variance  of  K  explained  by  t)ie  polynomials  (0.69)  is  much  smaller  than 

t'tiat  for  R  and  T.  :  (lowever,  the  calculation  of  is  r.uch  less  sensitive 
k  k  3 

to  W  than  to  R,  and  T,  and  therefore  it  is  felt  that  the  accuracy  in  es- 
k  k 

timating  K  from  the  bi-quadratic  polynomial  is  adequate.  After  determin¬ 
ing  appropriate  values  cf  W,  appropriate  values  of  R  and  T  are  found 

K  K 

from  Eqs.  (12)  and  (13). 

With  diffuse  radiation  .conditions,  the  appropriate  basic  state  val¬ 
ues  of  R  and  T,  are  not  functions  of  cos  z  and  therefore  are  not  deter- 
k  k 

mined  from  the  cubic  polynomial.  The  diffuse  basic  state  values  are 

listed  in  Tables  2  and  3  under  column  D  for  layers  2  and  3.  In  ..sing 

values  of  D  in  a  layer  containing  fractional  cloudiness,  we  assume  that 

W  =  1.  Therefore,  R^  and  T^  for  this  layer  are  the  proportionally 

weighted  averages  of  the  clear  and  overcast  valuec  of  D. 

After  obtaining  the  seven  parameters  (R^,  R^,  R^,  ,  T^ ,  T^,  T^) 

Eq .  ^9)  is  solved  for  X  ,  which  is  expressed  in  terms  of  X  ,  the  extra- 

J  O 

terrestrial  radiation  per  unit  horizontal  surface.  To  express  in 
terrris  cf  energy  received  at  t)ie  ground,  we  must  ev.iluate  X^. 


I  i 

I  - 
I  A 

i  -| 


"o  =  ^  ‘d’' 


(22) 


where  S  is  the  solar  constant,  d  the  distance  and  d  the  mean  distance  of 
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66 

the  earth  from  the  sun.  According  to  Hoyt  (1981)  ,  the  best  mean  value 

2  -  2 

tor  cue  qviiet  sun  solar  constant  is  1369.2  W/m  .  (ci/d)  is  a  function  of 
the  ellipticity  of  the  earth's  orbit  and  the  position  of  the  earth  in  its 
orbit  around  the  sun.  It  may  be  expressed  as 

(d/d)  =  [  1.000140  +  0.016726  cos  27T(JD-2)  I  (23) 

365 

■.  e  JD,  the  Julian  day,  is  1  on  January  1  and  365  on  December  31. 


CUBIC  POLYNOMIAL  COLFFICIENTS 

layer)  and  (overcast  layer) 

FOK  REFLE. 

,  k  =  1, 

CTIVITV  (R^) 

2,  3. 

:  r.  (clear 
k 

=  a  +  a,  cos 
0  1 

2 

z  +  cos  z 

+  cos 

3 

c 

a 

a , 

d 

0 

1 

2 

3 

.12395 

. 34765 

.  39478 

-.14627 

"2 

.15325 

. 39620 

. 42095 

- . 14200 

>^3 

.15946 

.42185 

. 48800 

-.18493 

"3 

(F/K) 

.27436 

.43132 

.  26920 

-.00447 

Pi 

(THIN  Ci/Cs) 

.25674 

.18077 

-.21961 

. 25272 

Pi 

(THICK  Ci/Cs) 

.60540 

.55142 

-.23389 

.43648 

P2 

(As/Ac) 

.66152 

. 14863 

-.08193 

. 13442 

P3 

(Sc/St) 

.71214 

.15033 

.00696 

. 03904 

P3 

(Cu/Cb) 

.67072 

.13805 

-.10895 

. 09460 

NOTES;  (F/K)  is  used  for  the  clear  layer  reflectivity  of  the  bottom 

layer  when  foq  and/or  smoke  are  reported. 

High  clouds  (layer  1)  -  any  form  of  cirrus  or  cirrostratus  or 
cirrocumulus  -  are  distinguished  only  by  thin  or  thick. 

Middle  clouds  (layer  2)  -  all  forms  of  middle  clouds  are  governed 
hy 

Low  clouds  (layer  3)  -  we  distinguisli  only  between  stratiform  clouds 
( stratocumulus ,  stratus,  nimbostratus)  and  convective  clouds  (cumu- 
Jus,  cumulonimbus). 


66 


Hoy  t , 


u. 


V.  , 


1981  : 


Personal  communication. 
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TABLE  11. 


CUBIC  POLYNOMIAL  COEFFICIENTS  FOR  TRANSMISSIVITY  (T  ) 


(clear 


layer)  AND  (overcast  layer)  ,  k  = 

1,  2,  3. 

T^  =  bo  +  b^  COB 

^  2 
z  cos 

Z  b3 

3 

cos  z 

b 

b. 

b. 

b 

o 

1 

2 

3 

.76977 

.49407 

-.44647 

.11558 

"2 

.69318 

.68227 

-.64289 

.17910 

S 

.68679 

.71012 

-.7x463 

.22339 

(F/K) 

.55336 

.61511 

-.29816 

-.06663 

(THIN  Ci/Cs) 

.63547 

.35229 

.08709 

-.22902 

(THICK  Ci/Cs) 

.264(8 

.66829 

.24228 

-.49357 

^2  (As/Ac) 

.19085 

.32817 

-.08613 

-.08197 

(Sc/St) 

.13610 

. 29964 

-.14041 

.00952 

I  (Cu/Cb) 

.17960 

.34855 

-.14875 

.01962 

NOTES:  tj(F/K)  is  used  for  the  clear-layer  transmissivity  of  layer  3  when 
foq  and/or  smoke  are  reported. 

Remaining  comments  same  as  those  in  Table  10. 


TABLE 

12.  BI -(QUADRATIC 
CLOUD  TYPE  AS 

POLYNOMIAL 

A  FUNCTION 

COEFFICIENTS  FOR  THE 
OF  CLOUD  FRACTION  (f 

WEIGHT  (W) 

:  )  AND  COS 
k 

FOR  EACH 

Z. 

W  -  C  +  C,  cos  z  +  f 

O  1  2 

k  ^k 

cos  z  +  C. 

4 

2 

cos  z  +  C 

c 

C, 

C, 

C. 

C.. 

o 

1 

2 

3 

4 

5 

THIN 

Ci/Cs  0,675 

-3,432 

1.929 

0.842 

2.693 

-1.354 

THICK 

Ci/Cs  1.552 

-1.957 

-1.762 

2.067 

0.448 

0.932 

As/Ac 

1.429 

-1.207 

-2.008 

0.853 

0.324 

1.582 

Sc/St 

0.858 

-1.075 

-0.536 

0.750 

0.322 

0.501 

Cu/Cb 

2.165 

-1.277 

-3.785 

2.089 

-0.387 

2.342 
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6.  ANALYSIS  OF  ERROR  AND  CONCLUSIONS 


Error  analyses  are  derived  separately  for  the  13  development  stations 
and  the  12  teat  stations  as  well  as  for  all  stations  combined.  Figure  3 
shows  the  geographical  distribution  of  the  stations  used  for  development 
of  the  reflectivity,  absorptivity  and  transmissivity  coefficients  as  well 
as  the  stations  reserved  for  testing.  Although  the  SOLMET  tapes  contained 
data  for  Omaha,  NE  and  Stephenville ,  TX,  these  stations  were  not  used. 

The  cloud  information  for  Omaha  was  missing  and  the  data  sample  for 
Stephenville  (which  is  close  to  Fort  Worth)  was  too  small  to  be  useful. 

Although  the  stations  are  designated  as  development  and  test  sta¬ 
tions  only  the  so-called  uniform  cloud  states  were  used  in  developing 


the  coefficients.  These  states  included  clear  skies  (that  is,  less  than 
0.1  cloud  cover)  and  overcast  conditions  where  the  first  cloud  layer 
visible  from  the  ground  was  overcast.  Of  the  total  number  of  complete 
hourly  observations  available  at  the  development  stations  (660S04) ,  less 
than  35  percent  consisted  of  uniform  cloud  states  as  defined  above.  The 
remainder  consisted  of  multi-layer  cloud  states  and  fractional  cloud 
states.  From  the  latter  group,  single-layer  fractional  cloud  states  were 
used  to  determine  the  weighting  coefficients.  These  data  represent  less 
than  20  percent  of  the  total  development  sample.  Thus  in  testing  the  mod¬ 
el  within  the  development  sample,  more  than  65  percent  of  the  data  is 
partially  independent  and  nearly  50  percent  of  the  data  is  wholly  inde¬ 
pendent  of  the  data  used  to  develop  the  computational  details.  Neverthe¬ 
less,  we  will  refer  to  this  sample  as  the  developmental  sample.  The  test 
sample,  consisting  of  a  comparable  number  of  hourly  observations  (579839), 
is  of  course,  completely  independent. 

Figures  5  and  6  show  res[jcctively ,  the  root  mean  sguure  error  (RMSF) 
and  the  bias  in  the  fraction  of  the  extra-terrestrial  solar  radiation 
transmitted  to  t)io  tjround  at  uacti  of  the  stations.  kMSE  is  defined  as 


(RI-1.SF)  .  = 


"  ,  N. 

=  L-  (X  -  X 


wliere  N  is  the  total  iiumlier  of  tiout  ly  obsc  rval  i  ons  af  '-.laUion  X  is 

1  n 

the  computed  value  oi  the  total  insolation  r''coivcd  .jt  the  ground  and  X 


around 


TEST  STATIONS 


is  the  observed  value.  Both  X  and  X  are  expressed  as  fractions  of  X„, 

n  n  ‘  0' 

the  solar  radiation  outside  the  atmosphere  per  unit  horizontal  surface. 


The  bias  is  simply 


N  . 

.  -  ’ 

=  i-l 


(X  -  X  ) 
n  n 


The  combined  values  of  RMSE  for  the  development  stations,  test  stations  and 
all  stations  shown  in  Part  A  of  Table  13  are  obtained  by  combining  the  sums 
of  the  squares  of  all  of  the  individual  errors.  That  is, 


r-  N . 

=  i  y-  y  (X  -  X 

N  n  n  . 

3=p  1-1  1 


where  p,  q,  and  N  are  respectively  1,  13,  660504  for  the  development  sta¬ 
tions-,  14,  25,  579339  for  the  test  stations;  and  1,  26,  1240343  for  all 
stations. 

The  slight  improvement  in  both  RMSE  and  BIAS  (see  Part  A  of  Tables  13 
and  l<i)  for  the  test,  stations  as  compared  with  the  development  stations, 
iilthooqh  ui'if  xiioote'i ,  iiinioao':edly  represents  a  sampling  iiuttuation  due  to 
t/ic  different  mix  of  otationt-  in  the  two  samples.  However,  xt  is  grati¬ 
fying  to  find  the  teui  res-jitc  no  worse  than  the  partially  dependent  sample. 

Aithouq;'  there  is  an  overall  positive  bias  of  about  1  percent,  it  does 
not  apjiaar  wortliwr.ile  to  modify  the  coefficients  inasmuch  as  the  bias  for 
ti'ie  test  stations  is  essentially  zero.  Furthermore,  no  geographical  pat- 
ter\i  of  bias  can  be  disoernod  from  Figure  6, 

RMClt;  ana  FIAo  were,  evaluated  as  functions  of  fractional  cloud  cover, 
nolaj'  zenith  angle,  the  number  and  distribution  of  cloud  layers  present 
and  ueathei  oanditit^ns.  The  results  are  shov/n  in  Tables  13  and  14.  N.ot 
uneHpsct.odly ,  the  smallest  error  occurs  with  clear  skic.s.  The  error  in¬ 
creases  v/itJ,  increas  ing  cloud  fraction  up  to  ninc/tenths  cloud  coverage 
afid  js  s-onewhat.  smaller  for  overcast  shies.  Part  of  this  distribution  of 
r^MSE  is  due  to  the  fact  that  the  clear  and  overcast  categories  are  essen¬ 
tially  without  ohaetver  eiroz'  in  estimating  sky  cover.  Judging  from  tlie 
di  .str  ibutioti  of  error  ,  this  component  of  error  appears  to  account  for 
about  20  percent  of  the  V.otol  RMSi:.  The  remaining  sources  of  error,  which 
will  bn  discua.sec!  later  in  this  section,  are  obviously  increasingly  prom- 
irient;  with  j.ncroa'dr.a  sky  envoi. 
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TABLE  14.  AVERAGE  BIAS  OF  MODEL  CALCULATION  OF  .  EMISSION  OF 

SOLAR  RADIATION  TO  GROUND  FOR  VARIOUS 


A.  FRACTIONAL  CLOUD  COVER 


o 

o 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

TOTAL 

DEV 

-.003 

.005 

.004 

.004 

.011 

.018 

.  022 

.027 

.039 

.057 

.028 

.018 

TEST 

-.001 

-.006 

-.009 

-.011 

-.011 

-.010 

-.006 

-.001 

.008 

.022 

.005 

.000 

ALL 

-.002 

-.000 

-.001 

-.002 

,003 

.007 

.010 

.015 

.  026 

,042 

.018 

.010 

B. 

COS  Z 

(INTERVALS  CENTERED 

ON) 

0.05 

0.15 

0.25 

0.  35 

0.45 

0.55 

0.65 

0.75 

0.85 

0.95 

DEV 

.018 

.027 

.011 

.  004 

.011 

.014 

.028 

.027 

.020 

.029 

TEST 

.013 

.012 

-.006 

-.014 

-.007 

-.006 

-.008 

-.007 

.002 

.004 

ALL 

.016 

.020 

.003 

-.004 

.002 

.004 

.019 

.018 

.011 

.019 

C.  lAYERED  CLOUD  DISTRIBUTION 


L 

M 

H 

LM 

LH 

MH 

LMll 

CLEAR 

DEV 

.011 

.015 

.000 

.041 

.036 

.028 

.070 

-.003 

TEST 

-.016 

-.004 

-.017 

.018 

.012 

.007 

.043 

-.001 

ALL 

.001 

.006 

-.003 

,032 

.027 

.017 

.058 

-.002 

D .  WEATHER 

NO  WEATHER 

PRECIPITATION 

FOG /SMOKE 

DEV 

.017 

.011 

.014 

TEST 

.001 

-.021 

-.012 

ALL 

.010 

-.001 

.001 

1  To  place  these  RMSE's  in  perspective  we  can  compare  them  with  the  ob¬ 

served  standard  deviation  (n)  of  for  comparable  categories.  The  only 
strictly  comparable  category  for  which  O  is  available  is  that  for  clear  skies. 
The  O  for  this  category  is  0.078,  which  is  to  be  compared  with  the  RMSL  of 
0.083.  This  means  that  for  this  category  the  error  in  the  model  calculatit^n 
is  about  the  same  size  as  the  variability  of  the  observations.  Inasmuch  as 
the  model  yields  a  fixed  value  for  a  particular  set  of  conditions,  such  as 
clear  skies,  a  represents  a  lower  bound  for  RMSE.  The  small  difference  be- 
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tween  0.078  and  0.083  is  largely  accounted  for  by  the  small  bias  in  the 
calculation  for  this  category. 

Another  category  that  can  be  compared,  although  not  with  the  same  pre¬ 
cision,  is  that  for  overcast  skies.  In  this  case  RMSE  is  0.146  whereas  the 
average  standard  deviation  of  the  observations  for  overcast  skies  is  0.129. 
However,  the  RMSE  applies  to  all  overcast  sky  states,  including  states  with 
multiple  cloud  layers,  whereas  the  O  pertains  only  to  overcast  sky  states 
where  the  first  cloud  layer  visible  from  the  ground  is  overcast.  IE  we  ex¬ 
amine  Part  C  of  Table  13  we  find  that  RMSE  is  larger  for  conditions  with 
multiple  cloud  layers.  In  view  of  this  variation  of  RMSE  with  multiple 
cloud  layers,  the  small  difference  between  the  RMSE  of  0.146  and  the  O  of 
0.129  would  undoubtedly  be  diminished  further  if  the  RMSE  category  were 
available  for  the  same  "uniform"  or  single-layer  overcast  conditions  for 
which  the  0  applies.  In  addition,  we  shall  show  that  the  "true"  RMSE  is 
even  smaller  than  is  indicated  in  Table  13. 

From  Part  B  of  Table  13  it  is  apparent  that  except  for  the  cases  where 
the  sun  is  close  to  the  horizon,  RMSE  varies  little  with  cos  z.  This  re¬ 
sult  indicates  that  the  variability  of  X  with  cos  z  is  being  properly 

n 

handled  by  the  model.  The  larger  error  for  the  lov.-est  cos  z  category  is  a 
consequence  of  the  manner  in  which  the  SOIJ4ET  data  were  recorded  and  does 
not  necessarily  indicate  poorer  model  performance. 

As  already  indicated,  RMSE  is  larger  for  multi-layered  cloud  states. 
This  result  is  not  unexiiected  inasmuch  a.s  with  more  complicated  cloud 
states  there  is  more  opportunity  for  observer  error.  The  last  part  of 
Table  13  shows  the  variation  of  RM8E  with  three  weather  states.  While  the 
error  is  larger  when  there  is  precipitation  as  com()ared  with  no  weather  or 
smoke  or  fog,  the  difference  is  not  very  largo.  Furthermore,  the  RMSE  for 
firecipitation  is  almost  the  same  as  that  for  overcast  skies.  Thus  it  aji- 
pears  that  the  presence  of  precitiitation  does  not  in  itself  appreciably  in¬ 
crease  the  RMSE  above  that  which  is  expected  for  overcast  skies  witliouL 
precipitation. 

Table  14  shows  the  bias  as  a  function  of  the  .some  parameters  as  in 
Table  13.  The  overall  bias  is  about  1  percent  positive  indicating  that  the 
calculated  radiation  reaching  rhe  ground  is  sliqtitly  larger  than  the  ob¬ 
served,  The  errors  appear  to  be  larger  for  the  so-called  developmental 
sample  than  for  the  test  samiilc,  but  this  difference  is  undoubtedly  for- 
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tuitous.  There  appears  to  be  an  increase  in  bias  with  increasing  fractional 
cloud  cover.  However,  the  bias  for  overcast  skies  is  less  than  that  for 
both  0.8  and  0.9  cloud  cover.  Since  the  calculated  insolation  for  frac¬ 
tional  cloud  cover  is  a  function  of  that  for  clear  and  overcast  states,  it 
seems  likely  that  the  apparent  relationship  between  RMSE  and  cloud  fraction 
is  fortuitous.  tt  is  possible  that  improper  values  of  the  weighting  factor 
W  are  contributing  to  the  variation  of  RMSE  with  cloud  cover,  but  if  the 
W  values  were  an  important  contributor  it  might  be  expected  that  the  ap¬ 
parent  relationship  of  RMSE  and  cloud  cover  would  be  enhanced  in  the  test 
sample.  It  is  obvious  that  such  is  not  the  case.  In  any  event,  the  bias 
in  the  test  sample  it  not  in  the  overall  seunple  is  small  enough  to  be  ig¬ 
nored  . 

From  Part  B  of  Table  14  it  is  apparent  that  there  is  no  consistent 
variation  of  bias  with  zenith  angle.  Therefore,  here  too  we  may  assume 
that  the  effects  of  zenith  angle  are  being  handled  appropriately. 

From  Part  C,  it  appears  that  it  might  be  possible  to  decrease  the  bias 
by  modifying  the  coefficients  for  multi-cloud  layers,  particularly  when 
clouds  are  present  in  all  three  layers.  it  is  not  obvious,  however,  that  the 
benefits  to  be  derived  would  be  worth  the  added  complexity  in  the  computer 
code,  especially  since  Part  D  indicates  that  the  bias  does  not  seem  to  de¬ 
pend  heavily  upon  the  state  of  present  weather. 

Judging  from  the  results  shown  in  Tables  13  and  14,  along  with  the 
limited  comparisons  with  comparable  standard  deviations,  both  the  coeffi¬ 
cients  derived  from  the  developmental  sample  and  the  procedure  for  calcu¬ 
lating  the  solar  radiation  received  at  the  ground  produced  successful  re¬ 
sults  on  completely  independent  data:  that  is,  there  was  no  loss  in  ac¬ 
curacy  with  independent  data.  This  is  not  to  say  that  both  the  coefficients 
and  the  methodology  cannot  be  improved.  However,  questions  on  the  desir¬ 
ability  of  such  improvements  in  terms  of  their  costs  are  beyond  the  scope 
of  this  report.  Nevertheless,  it  is  useful  to  examine  the  sources  of  er¬ 
ror  and  to  indicate  i>ossible  avenues  of  improvement.  Furthermore,  as  in¬ 
dicated  above,  we  shall  show  that  the  "true"  error  (and  for  that  matter 
the  "true"  standard  deviation)  is  less  than  that  indicated  in  Table  13. 

The  {.irincipal  sources  of  error  in  RMSE  can  be  classified  as  follows: 

1.  Observer  error 
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2.  Errors  due  to  lack  of  information 

3.  Errors  due  to  lack  of  representativeness 

4.  Measurement  error 

5.  Errors  due  to  model  simplicity 

We  have  already  alluded  to  observer  error.  This  arises  largely  from 
the  human  observer's  inability  to  integrate  accurately  the  fractional  cloud 
cover  of  an  individual  cloud  layer  except  for  those  cases  where  the  layer 
is  overcast  or  the  sky  is  clear.  In  addition,  observer  error  includes 
cases  where  the  cloud  type  is  inappropriately  identified.  There  is  little 
that  can  be  done  with  the  solar  insolation  model  to  reduce  this  source  of 
error  as  long  as  all  of  the  cloud  and  weather  information  is  derived  from 
standard  ground-based  observations. 

Another  important  source  of  error  due  to  the  ground-based  nature  of 
the  meteorological  observations  is  the  lack  of  cloud  information  for  layers 
above  the  lowest  overcast  layer.  A  possible  means  of  reducing  this  uncer¬ 
tainty  and  thereroro  minimizing  this  source  of  error  would  be  to  make  use  of 
available  cloud  analyses  (such  as  the  3-D  Nephanalysis  of  the  Air  Weather 
Service  ulobal  weather  Central)  which  incorporate  important  additional 
sources  of  information  such  as  satellite  information.  However,  one  of  the 
principal  virtues  of  this  solar  insolation  model  is  that  it  depends  only 
upon  routine  ground-based  observations. 

One  source  of  error  which  is  present  in  the  results  shown  in  Table  13, 
but  which  presumably  would  not  be  present  in  practice,  concerns  the  ground 
albedo.  Lacking  specific  information  on  the  ground  albedo  within  the 
SOLMET  data  sample,  a  fixed  climatological  mean  value  of  0.15  was  used  for 
all  stations  and  all  seasons.  in  practice,  however,  it  can  be  expected 
that  there  would  be  more  information  available  on  ground  albedo.  The  use 
of  such  information  could  only  diminish  RMSE,  but  of  course  would  have  no 
effect  on  O.  In  those  situations  in  which  the  actual  ground  albedo  dif¬ 
fered  significantly  from  0.15,  the  computed  solar  radiation  at  the  ground 
could  be  appreciably  modified  by  the  use  of  the  appropriate  ground  albedo. 
For  example,  with  a  low  overcast  and  high  sun  elevation  is  about  0.30 

with  R  =  .15  but  X  increases  to  0.45  X..^  with  R  =  .65.  For  the  same  con- 
g  3  0  ^  g 

ditions  but  with  low  solar  elevation  angle,  is  less  than  0.12  with 

R  =  .15  but  greater  than  0.19  X  with  R  ^  .65. 
q  0  g 
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Another  source  of  error  in  the  present  results  which  would  not  be 
present  in  practice  concerns  the  method  of  data  tabulation  in  the  SOLMET 
tapes.  Some  of  the  data,  such  as  the  pyrheliometric  data,  are  integrated 
over  an  hour,  but  other  data,  such  as  the  standard  meteorological  data, 
refer  to  a  specific  time  period.  In  situations  where  the  weather  is  chang¬ 
ing  rapidly,  the  weather  observations  for  the  assigned  hour  may  not  be  rep¬ 
resentative  of  the  insolation  measurement  for  that  hour. 

A  major  source  of  error  concerns  the  basic  solar  insolation  data.  Al¬ 
though  gross  errors  were  presumably  eliminated  from  the  SOIMET  tapes,  many 
inconsistencies  and  anomalies  remain  in  the  data. 

All  of  the  above  sources  of  error  contribute  to  the  indicated  RMSE  in 
Table  13,  but  are  model  errors;  that  is,  they  are  not  true  errors  as 

far  as  the  model  is  concerned.  How  much  of  the  indicated  RMSE  is  due  to 
such  errors  is  impossible  to  say,  but  they  must  be  appreciable.  Many  of 
these  sources  of  error  also  contribute  to  the  observed  standard  deviation 
so  that  it  is  clear  that  the  "true"  standard  deviation  is  also  smaller 
than  is  indicated  by  the  data. 

There  are  riodel  errors,  of  course,  and  these  are  largely  due  to  misrep¬ 
resentations  and  simplifications  of  the  radiation  physics  by  the  modeling 
assumptions.  Prominent  among  these  simplifications  are  (1)  the  truncation 
errors  introduced  by  representing  a  continuous  atmospheric  medium  with 
only  three  discrete  layers,  (2)  the  simple  and  arbitrary  treatment  of  dir¬ 
ect  and  diffuse  radiation,  (3)  the  assumption  that  the  (R,  T,  A)  coeffi¬ 
cients  are  the  same  for  upward  as  for  downward  directed  radiation,  (4)  the 
neglect  of  moisture  and  aerosols  as  a  function  of  time  and  place,  (5)  the 
implicit  neglect  of  seasonal  and  geog,. aphical  variations  by  using  the  same 
coefficients  for  all  stations  for  all  seasons,  (6)  the  use  of  a  limited 
number  of  basic  cloud  type  categories,  and  (7)  the  treatment  of  solar  rad¬ 
iation  foi:  all  practical  purposes  as  monochromatic. 

All  of  the  above  simplifications  and  approximations  can  be  reduced 
but  only  at  the  expanse  of  increased  model  complexity. 

In  view  of  the  errors  in  the  observational  data  needed  to  evaluate 
improved  (R,  T,  A)  coefficients,  such  model  improvements  do  not  appear  to 
be  warranted  at  the  present  time.  However,  it  is  also  appropriate  to 
point  out  that  in  spite  of  the  simplifications  and  approximations  the 
model  appears  to  estimate  the  solar  radiation  received  at  the  ground  with 
relatively  little  error. 
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